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ABSTRACT 


Two  theorems  were  derived.  First,  the  Initial  Response  Theorem 
describes  the  necessary  and  sufficient  conditions  for  a  series  compensator 
for  a  feedback  control  system  to  simultaneously  stabilize  the  system, 
and  cause  the  initial  system  response  to  achieve  prescribed  constraints. 

A  single  input-single  output,  continuous  time  linear  system  is  considered, 
with  a  delta  function  driving  any  rational  transfer  function  in  the 
s-domain  being  the  system  input.  Design  constraints  can  be  placed  on 
the  initial  response  value,  and  on  any  of  the  derivatives  (from  the  right) 
of  the  initial  response.  Second,  the  Initial  Response  Parameterization 
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gives  a  parameterization  of  the  complete  set  of  compensators  that  will 
meet  the  given  constraints  whenever  the  conditions  of  the  Initial 
Response  Theorem  are  met.  The  area  of  Youla,  Bongiorno  and  Jabr  (YBJ) 
Control  Theory  was  used  for  the  derivations,  but  first  required  a 
system  transformation  to  convert  the  initial  value  problem  to  one  of 
stability. 
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CHAPTER  I 


INTRODUCTION 


1.1  Objective 

The  results  of  this  thesis  are  two  theorems.  The  first  is  named 

/ 

the  Initial  Response  Theorem  (IRT),  which  describes  the  necessary  and 
sufficient  (N  &  S)  conditions  for  a  series  compensator  for  a  feedback 
control  system  to  simultaneously  stabilize  the  system,  and  cause  the 
initial  system  response  to  achieve  prescribed  constraints.  A  single 
input-single  output,  continuous  time  linear  system  is  considered.  The 
system  can  have  as  its  input  a  delta  function  driving  any  rational 
transfer  function  in  the  s-domain.  Design  constraints  can  be  placed 
on  the  initial  response  value,  and  on  any  of  the  derivatives  of  the 
initial  response,  i.e.,  as  time  approaches  0^  (0  from  the  positive 
side).  The  second  theorem  is  termed  the  Initial  Response  Parameterization 
(IRP),  and  it  gives  a  parameterization  of  the  complete  set  of  compensators 
that  will  meet  the  given  constraints  whenever  the  conditions  of  the 
IRT  are  met. 

The  motivation  for  this  work  was  the  intuition  that  the  shape  of 
the  transient  response  can  be  controlled  by  controlling  the  initial 
response.  Thus,  if  the  initial  response  and  the  asymptotic  stability 
can  be  simultaneously  controlled,  the  designer  will  have  a  tool  that 
can  help  meet  both  transient  response  requirements  and  stability 
requirements  simultaneously.  Examples  of  this  technique  are  given  in 
Chapter  IV. 


The  theorems  are  derived  in  essentially  two  stages.  First,  the 


given  system  is  modified  to  another  system.  Any  compensator  causing 
the  modified  system  to  be  stable  will  also  cause  the  original  system 
to  meet  the  desired  constraints.  Then,  the  relatively  new  area  of 
YBJ  control  theory,  which  is  introduced  in  Section  1.3,  is  applied  to 
the  modified  system.  This  leads  to  the  N  &  S  conditions,  and  the 
parameterization  describing  the  complete  set  of  compensators  that 
will  stabilize  the  modified  system,  and  therefore  also  meet  the 
constraints  given  for  the  original  system.  The  constraint  of  system 
stability  is  achieved  by  application  of  previous  results  in  YBJ  (10). 
the  initial  condition  constraints  are  met  via  the  application  of  new 
results  in  YBJ  derived  in  this  paper.  All  of  the  results  assume  that 
the  plant  considered  is  a  proper  rational  function. 


Figure  1.  Example  System 


As  an  example,  consider  the  system  shown  in  Figure  1.  The  input 
is  an  impulse,  which  is  represented  by  a  delta  function  driving  a 
transfer  function  of  value  1.  The  plant  has  a  transfer  function 


and  is  unstable.  It  is  desired  that  the  system  H  (s)  be  stable, 

Vl 

and  that  the  initial  response  V2(0^)  equal  3. 


The  parameterization  of  the  complete  set  of  compensators  that 
stabilize  the  feedback  loop  has  been  found  before  this  thesis  (10). 
Without  going  into  further  detail  at  this  point,  it  will  suffice  to 
state  that  this  parameterization  leads  to  c(s)  of  the  form 


c(s)  = 


<1.2> 


where  w(s)  is  any  stable  function. 


The  compensator  must  also  meet  the  constraint  V2(0^)  equal  3. 


The  system  has  the  transfer  function 

_  1  c(s)p(s)  <1.3> 

H(s)  -  1  i-irtflTpIiT- 

If  <1.1>  and  <1.2>  are  substituted  into  Equation  1.3,  then  after  some 


manipulation 


H(s)  =  -w(s)  ^^2  ill 


<1.4> 


(s+1) 


Also,  according  to  the  Initial  Value  Theorem 


lim  {sH(s)}  =  lim  {h(t)} 
s-^o"  t-K)+ 


<1.5> 


where  H(s)  is  the  Laplace  Transform  of  h(t).  Thus,  it  is  also  required 


lim  {-w(s)  =  3 

s^  (s+1)^ 


<1.6> 


in  order  to  meet  the  initial  response  constraint.  One  could  then  search 


"brute  force"  for  stable  w(s)  which  satisfy  Equation  1.6  in  order  to 
solve  the  problem,  and  could  achieve  particular  solutions  such  as 


<1.7> 


However,  the  IRP  in  Chapter  III  yields  the  complete  set  of  compensators. 
The  calculation  is  not  difficult  for  this  problem,  and  will  be  omitted 
here  to  keep  the  example  brief.  But,  when  applied,  the  IRP  yields 


w(s)  » 


(s+1) 


s+1 


<1.8> 


as  the  complete  parameterization,  where  e(s)  is  arbitrarily  stable,  and 
c(s)  is  again  described  in  Equation  1.2. 


1.3  YBJ  Control  Theory 

The  foundation  of  YBJ  control  theory  is  two  papers  published  in 
1976  by  Youla,  Bongiorno,  and  Jabr  (16,  17).  A  result  of  their  work 
was  the  complete  parameterization  of  the  set  of  stablizing  compensators 
for  a  multivariate  feedback  system,  based  on  a  new  approach  in  feedback 
system  design.  Since  then,  this  new  approach  has  led  to  similar 
parameter! zations  of  other  problems,  such  as  the  tracking  and  disturbance 
rejection  problems  (9,  10,  11),  and  also  to  results  in  optimization 
theory  based  on  the  parameter! zations  found  (16,  17). 

The  basic  approach  to  discovering  the  parameter! zations  (other 
than  the  original  stability  parameterization)  has  been  to  arrange  the 
mathematics  so  that  they,  too,  were  a  stability  problem.  Thus  similar, 
though  different  approaches  to  the  original  stabilization  problem  could 
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be  used.  This  then  is  the  reason  for  the  system  transformation  used 
in  the  derivation  of  the  theorems  in  this  paper,  as  will  be  seen  in 
the  following  chapters. 


The  key  to  the  YBJ  approach  is  to  use  a  "stable  fractional 
representation"  for  the  transfer  functions  considered.  This  is 
opposed  to  the  "polynomial  fractional  representation"  used  classically. 
For  example,  the  function  r(s)  would  classically  be  written  in  the  form 

r(s)  =  ^  <1.9 

where  p(s)  and  q(s)  are  polynomials.  But,  YBJ  requires  that  r(s)  be 
wri tten 

n.(s) 

where  n^(s)  and  d^(s)  are  both  stable,  and  have  no  common  (closed)  right 
half-plane  (RHP)  zeros,  including  infinity.  This  is  termed  "RHP 
coprimeness",  or  simply  "coprimeness".  Thus,  if  m(s)  is  a  Hurwitz 
polynomial  of  order  equal  to  the  order  of  r(s),  then 


and 


are  possible  representations  for  YBJ.  Then,  in  the  case  that  n^(s) 
and  d^(s)  are  RHP  coprime,  YBJ  dictates  that  there  must  exist  stable 
u^(s)  and  v^(s)  such  that 

u„(s)n„(s)  +  v^(s)d^(s)  =  1  <1.13> 


Figure  2.  General  System  Configuration 


Next  define 


n  (s) 

nt(s) 

y=)  =rT5r 

n 


<1.15> 


<1.16> 


where  n  (s)  is  stable  and  coprime  with  stable  d  (s),  and  where  n^(s)  is 
P  ^  ''n 

stable  and  coprime  with  stable  d.(s). 

^n 

With  the  above  definitions  in  mind,  the  following  theorom  is  a 
major  result  of  this  thesis. 

Initial  Response  Theorem 

For  the  system  in  Figure  2,  assume  that  r{p(s)}^0.  Also  assume 


The  input  t(s)  has  fewer  than  three  poles  at  s  =  0. 


If  the  above  condition  is  not  met, 

then  {(#poles  t(s)  at  s=0)-3}<{#zeros  p(s)  at  s=0}. 


<1.17> 


<1.18> 


Then,  a  set  of  compensators  c(s)  exists  that  will  simultaneously  stabilize 
H  (s)  (without  RHP  pole-zero  cancellations  between  p(s)  and  c(s)),  and 

meet  the  initial  condition  constraints  {v5(0''') ,. . .  ,v^(0''') }  =  {,/«,.  ...y,-} 


if  and  only  if  (iff) 

r{t(s)}  +  r{p(s)}  <  M+1 


<1.19> 


<1.20> 


d.(s)  and  d  (s)  are  RHP  coprime, 

j  ^ 

which  can  be  written  {d^{s),  d  (s))  =  1. 

z.  p 

EOT  (End  of  Theorem) 

The  proof  for  this  theorem,  as  well  as  the  resulting  parameterization 
for  the  complete  set  of  adequate  compensators  are  given  in  Chapter  III. 

The  conditions  given  in  equations  <1.19>  and  <1.20>  are  not  at  all 
obvious.  In  practice,  <1.20>  should  not  normally  pose  a  problem,  while 
<1.19>  will  often  be  the  more  restrictive  condition. 

Though  a  different  approach  will  be  used  in  Chapter  III  to  prove 
<1,19>,  some  intuition  can  be  given  (unrigorously)  for  its  necessity 
now.  It  can  be  shown  that  the  system  ,,(s)  has  a  relative  degree 

V2U1 

equal  to  or  greater  than  that  of  the  plant  p(s),  whenever  r{p(s)}  >  0. 
Therefore,  V2(s)  must  have  relative  degree 

r{V2(s)}  ^  r{t(s)}  +  r{p(s)}  <1.21> 

since  the  relative  degree  of  the  product  is  the  sum  of  the  relative 
degrees.  Also,  the  Initial  Value  Theorem  (IVT)  can  be  used  to  prove  that 
r{V2(s)}  =  M+1  <1.22> 

Thus,  equating  <1.21>  and  <1.22>  leads  to  <1.19>. 

The  assumptions  made  in  equations  1.17  and  1.18  are  used  to 
guarantee  that  a  compensator  meeting  the  n-th  derivative  constraint 
will  also  meet  the  0-th  through  {n-l)th  constraints.  However,  they  may 
not  be  necessary  conditions.  At  this  time,  a  set  of  necessary  and 
sufficient  conditions  have  not  been  found,  though  the  given  assumptions 
are  quite  unrestrictive.  If  they  do  pose  a  problem,  the  interested 


reader  is  referred  to  Chapter  III  for  further  details.  Such  a 
"guarantee"  will  be  useful  when  applying  the  parameterization  to  be 
derived  later. 

The  remainder  of  this  thesis  is  organized  as  follows.  Chapter  II 
develops  the  system  transformation,  and  develops  a  list  of  properties 
associated  with  the  transformation.  Chapter  III  applies  the  YBJ  theory 
to  the  transformed  system  to  discover  the  Initial  Response  Theorem  and 
the  resulting  parameterization.  Next,  Chapter  IV  contains  design 
examples  using  this  theory.  Finally,  Chapter  V  contains  a  list  of 
recommendations  for  further  research  on  this  topic. 


This  chapter  deals  with  the  development  of  a  system  transformation 
that  modifies  a  given  feedback  system  model  so  that  YBJ  control  theory 
can  be  applied  in  the  next  chapter.  The  transformation  allows  for 
simultaneous  system  stability  and  initial  response  constraints  to  be 
attacked  with  YBJ.  The  system  is  single  input-single  output  continuous 
time,  and  linear.  The  system  input  is  a  delta  function  driving  a 
rational  transfer  function. 

As  was  discussed  in  the  introduction,  the  approach  to  deriving 
parameterizations  for  compensators  using  YBJ  in  the  past  has  been  to 
look  at  each  problem  as  one  in  stability.  In  some  way,  the  algebra  of 
the  problem  was  arranged  so  that  stabilizing  some  quantity  would  produce 
the  desired  results.  Then,  the  YBJ  theory  could  be  used  to  find  the 
parameterization  of  compensators  that  would  indeed  stabilize  the  quantity 
in  question.  Thus,  the  intent  here  is  to  derive  a  general  system 
transformation  that  leads  to  particular  quantities,  which  if  stabilized, 
will  cause  the  constraints  to  be  met. 

The  problem  at  hand  may  have  many  constraints,  but  essentially 
only  two  types.  First,  the  feedback  system  must  be  stable.  Second, 
some  initial  value  constraints  must  be  met.  Therefore,  it  is  necessary 
to  derive  some  transformation  that  produces  two  quantities.  The  first 
will  stabilize  the  feedback  loop  when  stabilized,  and  the  second  will 
meet  a  desired  initial  condition  constraint  when  stabilized.  Thus,  the 
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basic  transformation  concepts  to  be  used  here  are  as  follows.  First, 
the  transformation  of  the  feedback  loop  preserves  the  properties  of 
stability,  so  that  stabilizing  compensators  in  the  transformed  domain 
also  stabilize  in  the  original  s-domain.  Because  of  this  property, 
the  stability  constraint  can  be  met  using  a  previously  developed 
stability  parameterization.  Second,  the  transformation  yields  a 
system  output  that  meets  one  of  the  initial  conditions  when  stabilized. 
Also,  the  transformation  must  have  an  inverse  transformation,  so  that 
the  answer  can  be  of  use  in  the  original  s-domain. 

2.2  Mathematical  Transformation 

The  first  step  in  developing  the  system  transformation  is  to 
develop  a  mathematical  transformation  that  maps  an  initial  derivative 
response  into  the  final  value  response  in  a  new  domain.  For  this 
transformation  let  t  and  s  represent  the  time  and  Laplace  domain 
variables  of  a  transfer  function  as  usual.  Then,  let  t  and  £  represent 
the  corresponding  time  and  Laplace  domain  variables  in  the  transformed 
domain.  As  such,  the  following  lemma  describes  the  mathematical 
transformation  that  performs  the  desired  mapping. 

Lemma  2.1 

Let  h'^(s)  be  the  n-th  derivative  of  a  transfer  function 
H(s)  =  L{h(t)}.  If  H^(s)  is  strictly  proper,  and  if  H^(l)  only  has 
poles  in  the  left  half  plane,  except  possibly  for  a  simple  pole  at 
s=0,  then 


I 


Tim  (£)}  =  lim  {h„(t)} 
s->0  "  t-x»  " 

=  lim  {h'^(t)}  =  lim  {sH^(s)} 

S-x» 

where  H„(s)  =  ^  H(i)  -  h’(0*) 


Proof 


The  n-th  derivative  of  a  transfer  function  H(s)  is  equal  to 
h'^(s)  =  s"h(s)  -  s"'^“’h’(0'^) 


n  >  0 


where  h^(0‘^)  =  lim  {L"^{H^(s)}} 


The  Initial  Value  Theorem  (IVT)  states  that  if  H{s)  is  strictly  proper 

lim  {h(t)}  =  lim  {sH(s)}  <2 

t-hO'*'  s-Ko 

Of  course  if  H(s)  is  not  strictly  proper,  no  finite  initial  value 

exists  since  there  is  an  impulse  at  t  =  0.  Then,  combining  <2.4>  and 

<2.6>  leads  to 

lim  {h'^(t)}  =  lim  {sH'^(s)} 
t^+  s-^ 

=  lim  {s"'^^H(s)  -  s""''h\0'*’)  <2 


if  h'^(s)  is  strictly  proper. 

Next,  define  the  transfer  function 


in  another  domain,  with  the  given  variables.  Then,  the  Final  Value 


Theorem  (FVT)  states  that  if  H^(s^)  has  no  (closed)  RHP  poles,  except 
for  possibly  pole  at  s  =  0  (see  (8),  p.  714),  then 


lim  {h  (t)}  =  lim  (s)} 
t-H»  "  s-^  " 


<2.9> 


At  this  point,  a  change  of  variables  is  necessary  to  force 
lim  {^^(t)}  =  lim  {h'^(t)} 


<2.10> 


This  link  can  be  accomplished  by  forcing  equality  in  the  Laplace  domain. 


by  forcing 


lim  {sH  (£)}  =  lim  {sH'^(s)} 
s-»0  s-K» 


=  lim  {s"'^^H(s)  -  s'’"''h^(0+)} 

S-+®  i  =0 


<2.11> 


If  ^  and  s  are  related  by 


<2.12> 


lim  i  =  0  =  lim  ^ 
s-Hx-  s->0 


<2.13> 


<2.14> 


Thus,  equating  the  values  inside  the  brackets  in  Equation  2.10  dictates 


<2.15> 


which  leads  to 


1  ui/r>  +  \ 


<2.16> 


Therefore,  invoking  Equations  2.12  and  2.16  force  the  equality  of 
Equation  2.11,  and  therefore  also  Equation  2.10,  if  the  restrictions 
described  for  the  IVT  and  the  FVT  are  met. 

EOP  (End  of  Proof) 


Lemma  2.1  then  is  a  mathematical  transformation  yielding  a  new 
transfer  function  whose  asymptotic  value  is  equal  to  the  initial 
n-th  derivative  of  the  original  transfer  function  response.  Later 
in  this  chapter  a  method  of  applying  this  to  the  feedback  system  shown 
in  Figure  2  will  be  developed  so  that  the  YBJ  theory  can  be  used  to 
meet  the  n-th  derivative  initial  condition  constraint.  Also,  Section 
2.5  will  prove  that  the  restrictions  on  the  IVT  and  the  FVT  will  always 
be  met  in  this  application. 


2.3  Inversion  Transform  and  Properties 


Before  applying  the  transformation  defined  in  Section  2.2  to  a 
feedback  system,  it  will  prove  useful  to  first  define  the  transformation 


T  and  its  inverse  T”"^  as  follows. 
H(s)  =  f  {H(s)}  =  H(i) 

{H(s)}  =  H  (i) 


<2.17> 


The  following  properties  associated  with  this  transformation  will  also 
be  useful.  The  proofs  for  these  properties  can  be  found  in  Appendix  A. 


Additive  Property 

H(s)  +  F(s)  =  T  {H(s)  +  F(s)} 

Note  that  this  property  also  implies  the  subtractive  property. 


<2.18> 


ncative  Pro 


H{s)-F{s)  =  T  {H{s)-F{s)} 

Note  that  this  also  implies  division 

Inverse  Transform  Prooert 


H(|-)  =  H(s) 

Transformation  Stability  Pro 


If  H(s)  is  stable  (or  unstable)  in  the  s-domain,  then  so  is  H{^)  in 
the  s-domain,  and  visa  versa. 

2.4  System  Transformation 

A  method  of  using  the  results  of  Lemma  1  to  transform  the  system 

in  Figure  2  (reproduced  in  Figure  3,  part  a)  will  be  derived.  This 

will  result  in  the  n-th  derivative  of  the  initial  ouput  of  the  system 

lim  {v5(t)}  <2. 

t-^'^  ^ 

to  become  the  asymptotic  ouput  (time  approaches  infinity)  of  the 
transformed  system.  This  model  will  then  be  modified  so  that  the  YBJ 
stability  criterion  may  be  used  to  meet  the  condtion  in  Equation  2.25 


Def i ne 

the  transfer 

functions 

r'(s) 

i=0 

2(0+) 

r;(i) 

n-1  1 

1=0  £ 

Lemma 

2.2  applies. 

Lemma  2.2 


Given  the  system  in  Figure  3,  part  a,  and  the  functions  described 
in  Equations  2.26  and  2.27,  the  system  in  Figure  3,  part  b,  will  have 
an  ouput  ^3(1)  with  the  property 

lim  {e’^(t)}  =  lim  {v5(t)}  <2.28> 

t-x»  t^'^  ^ 

Proof 


The  transfer  function  of  the  system  in  Figure  3,  part  a  can  be 
calculated  to  be 

Then,  invoking  Lemma  2.1  would  require  a  system  with  the  transfer 
function 


KM  =  4r7  t(T)  vi(0*) 


=  in<^) 


<2.30> 


via  Equation  2.2,  and  invoking  the  multiplicative  and  additive  properties 
as  well  as  Equations  2.26  and  2.27. 


The  system  in  Figure  3,  part  b,  has  the  transfer  function 

iCjlPllJ—  _  r.(s) 


■3n®0n 


eo^Onl-'  "n  -  l+c(ijp(sj 


<2.31> 


Equations  2.30  and  2.31  are  identical,  revealing  that  Figure  3,  part  b, 
is  in  fact  consistent  with  Lemma  2.1. 
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YBJ  theory  requires  a  system  which,  when  stabilized,  meets  the 
required  constraints.  The  property  characterizing  any  stable  system 


is  that  asymptotically  the  response  approaches  zero.  Considering  this, 
if  a  different  transfer  function  can  be  given  that  asymptotically 
equals  the  response  of  the  system  in  Figure  3,  part  b,  then  the 
difference  between  these  would  be  asymptotically  zero.  Therefore,  a 
composite  system  equal  to  this  difference  would  be  stable. 


A  different  transfer  function  with  the  asymptotic  response  V2(0'*') 

is  the  step  function  with  weight  V2(0‘'’), 

n 
2 


w(s)  = 


v;(0+) 


<2.32> 


Then,  the  difference  leading  to  stability  would  be  the  transfer  function 

<2.33> 


^  e  (s)  =  H-  -  (s)  -  w(^) 
®4n®0n“  ®3n®0n- 


But,  an  equivalent  model  for  this  system  would  be  to  add  w(^)  to 


r^(s^),  creating  a  new  function,  follows, 


rj^(^)  =  w(£)  +  r^l^Cs) 


v"(0")  n-1  1  .  ^ 

=  — - +  L  — 

s 


5=0  ^2"” 


'  L  -2 


vi(0+) 


<2.34> 


i=0  s' 

This  then  leads  to  the  following  lemma. 


Lemma  2.3 

For  the  system  in  Figure  3,  part  a,  let  t(s),  p(s),  and  the 
initial  response  constraints  {yg,  y^, . . .  ,yj^,. . .  ,yj }  (y.j  is  finite,  real) 
be  given.  Then,  substitute  the  initial  response  constraints  into  >^^^(2), 


so  that 


r  (s)  =  E 


i=0  s 


n-i+1 


<2.35 


Then,  the  only  compensators  that  will  meet  the  constraint  V2(0'*’)  = 

which  may  possibly  also  meet  the  constraints  {V2(0'''),. . . ,  V2~^(0''')} 

=  {yn»---»  y„  i}*  those  which  stabilize  H-  -  (s)  in  Figure  3, 

^  ®4n®0n 

part  c,  where  is  written  per  Equation  2.35. 


Proof 


Transform  the  system  in  Figure  3,  part  a,  to  the  one  in  Figure  3, 


part  b.  Then,  in  accordance  with  Lemma  2.2 

lim  {e,(t)>  =  lim  {v5(t)}  =  v5(0‘’') 

^  t->0+  ^ 

Also,  let 

W(..  .  ^ 

lim  {w(t)}  =  y_ 


<2.36 


<2.37 


and  define  H-  -  (£)  by  Equation  2.33.  Then,  H-  -  (^)  will  be 
®4n®0n  ^40*^00 

stabilized  only  if  H—  -  (s)  {which  is  equal  to  e’,(s)}  and  w(s)  are 

®3n®0n“ 


asyptotically  equal.  This  of  course  requires  that 

y„  =  v5(o*) 


by  equating  Equations  2.36  and  2.37. 


<2.38 


Therefore,  consider  a  compensator  c(^)  which  leads  to  the  system 
in  Figure  3,  part  a,  with  inital  response  values  {v^lO'*') , . . . ,  V2"^(0'‘')} 
Then,  use  these  values,  and  let  V2(0+)  =  y^  for  Equation  2.34.  As 
such,  if  Hr  T  (s)  in  Figure  3,  part  c,  is  stable,  then  Equation  2.38 


Note  that  the  above  argument  does  not  consider  the  constraints 


^  compensator  may  or  may  not  meet  these 

constraints.  However,  if  it  does,  then  obviously  these  constraints 

can  be  substituted  for  {V2(0''') ,  v^’^CO'*’)}  in  Equation  2.34.  Thus, 

making  this  substitution,  and  invoking  Equation  2.38  into  Equation  2.34, 

leads  to  Equation  2.35.  Stabilizing  the  system  in  Figure  3,  part  c, 

using  Equation  2.35  therefore  leads  to  a  compensator  that  meets  the 

constraint  y^,  and  may  (no  guarantee  yet!)  meet  {yg,...,  Also, 

if  it  does  cause  the  system  to  meet  {yn»-«*»  y„  il*  then  H-  -  (s)  must 

^  ®4n®0n“ 

be  stable  when  defining  r^{^)  according  to  Equation  2.35. 
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Figure  4.  Final  System  Transformation 


For  the  system  in  Figure  4,  define 
tn(s)  =  {t^(s)}  =  s""'2t(s) 


r„(s)  =  r'  lr„(s))  = 


<2.39> 


per  the  transformation  in  Equation  2.17,  using  the  inverse  transform 


and  multiplicative  properties,  and  per  Equations  2.26  and  2.35.  Then 
the  following  theorom  applies. 

Theorem  2.1 

For  the  system  in  Figure  3,  part  a,  let  t(s),  p(s),  and  the 
initial  response  constraints  {yg,  y^,...  y^....,'yj}  be  given,  and 
define  t^(s)  and  >^p(s)  per  Equation  2.39.  Then,  the  only  compensators 
that  will  meet  the  constraints  v^CO'*’)  =  y^,  which  may  possibly  also  meet 
the  constraints  {v^CO"^),  v^CO"^),...,  =  {yg,  y^,... 

are  those  which  stabilize  H  (s)  in  Figure  4. 

®4n®0n 

Proof 

The  additive  and  multiplicative  properties  dictate  that  transforming 
a  system  from  the  s-domain  to  the  ^-domain,  and  visa  versa,  may  be 
accomplished  by  transforming  the  individual  transfer  functions  within 
the  system  separately.  Therefore,  the  system  in  Figure  4  is  the  s-domain 
equivalent  to  the  ^-domain  system  in  Figure  3,  part  c.  Also,  by  the 
stability  transformation  property,  these  systems  are  either  both  stable, 
or  both  not  stable.  Therefore,  since  the  statement  in  Lemma  2.3  relies 
only  on  the  stability  of  the  s^-domain  system,  an  equivalent  statement 
can  be  made  concerning  the  system  in  Figure  4. 

EOP 

With  Theorem  2.1  in  hand,  the  system  in  Figure  4  can  be  used  for 
the  application  of  the  YBJ  theory  to  find  the  stabilizing  compensators 
that  meet  the  initial  response  constraints.  This  will  be  accomplished 
in  Chapter  III. 
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2.5  Restrictions  Due  to  IVT  and  FVT 
The  system  in  Figure  4  is  the  model  from  which  a  complete  stabilizing 
parameterization  is  to  be  derived  in  the  next  chapter.  However,  this 
model  is  based  upon  the  mathematical  transformation  in  Section  2.2, 
which  incorporates  the  IVT  and  FVT.  At  that  time  it  was  mentioned 
that  both  of  these  theorems  have  restrictions  on  their  application.  The 
intention  here  is  to  prove  that  this  model  will  always  be  within  the 
bounds  of  these  restrictions,  so  that  a  complete  parameterization  can 
in  fact  be  derived  based  upon  it. 

First  consider  the  IVT,  which  states  that  if  H(s)  is  strictly 
proper,  then 

lim  {h(t)}  *  lim  {sH(s)}  <2.40> 

t'+O'*’  S-x” 

In  effect,  this  (strictly  proper)  says  that  there  must  be  a  finite 
initial  value  in  order  to  apply  the  theorem.  If  H(s)  is  not  strictly 
proper,  the  initial  value  is  undefined  (infinite),  and  therefore  equality 
cannot  hold,  even  though  the  IVT  will  also  yield  an  infinite  initial 
value.  Thus,  this  restriction  simply  requires  that  the  constraints 
{yg,  yj)  be  finite. 

The  FVT  states  that  if  H^(l)  has  poles  that  lie  entirely  within 
the  left  half  plane,  except  possibly  for  a  simple  pole  at  s=0  (see 
reference  (8),  p.  714),  then 
lim  {h  (0)  =  lim  {sH  (^)} 

t-Ko  "  " 


<2.41> 


For  the  system  model  here. 


<2.42> 


«.{!)  =  fV  i  (s)  -  r;(s) 

''2n®0n  " 

=  fV  ?  ^ 

''2n®0n  "  ^ 

~  ^  e  (s)  >  ^ 

®4n®0n  ^ 

and  —  (s)  will  always  be  stabilized.  Thus,  since  F_(£)  is  the 
®4n®0n  " 

sum  of  proper  functions,  it  will  be  proper.  Also,  its  poles  will  be 
those  of  its  components,  unless  some  are  cancelled.  Therefore,  it  will 
be  stable,  except  for  possibly  a  simple  pole  at  s=0,  so  the  FVT 
restrictions  will  always  be  met  also. 


2.6  Summary 

Stabilizing  the  system  model  in  Figure  4  has  been  shown  to  be 
N  &  S  to  guarantee  that  the  n-th  derivative  initial  condition  be  met, 
where  the  0-th  through  (n-l)th  initial  conditions  could  possibly  also 
simultaneously  be  met.  It  will  be  shown  in  the  next  chapter  that 
meeting  minor  restrictions  will  guarantee  the  first  n  constraints  are 
met  simultaneously  when  meeting  the  n-th  is  accomplished. 

Another  feature  of  this  model  is  that  the  feedback  loop  portion 
of  the  system  is  identical  to  that  in  the  original  system.  This  is 
important  since  it  is  also  desired  to  stabilize  this  loop,  and  convenient 
because  previous  work  in  YBJ  (10)  has  already  solved  this  problem. 


CHAPTER  III 


INITIAL  RESPONSE  THEOREM  AND  PARAMETERIZATION 


3.1  Introduction 

The  IRT  and  IRP  are  derived  in  this  chapter.  The  proceeding 
derivations  are  based  on  the  system  in  Figure  4  of  Chapter  II. 

3.2  Definitions  and  Properties 

Some  general  definitions  and  properties  from  reference  (10)  are 
needed  for  the  work  in  this  chapter. 

Definition  3,1 

Given  two  transfer  functions,  x(s)  and  y(s),  then  y(s)  divides  x(s)  if 


where  y'(s)  is  stable.  Notationally,  this  is  written  yjx. 

Definition  3.2 

Two  stable  transfer  functions  x(s)  and  y(s)  are  RHP  coprime,  or  simply 
coprime,  if  they  have  no  common  RHP  zeros.  Notationally,  this  is  written 
{x,y)  =  1. 

Definition  3.3 

The  transfer  functions  x(s)  is  said  to  be  mini  phase  if  it  is  stable, 
and  has  a  stable  inverse. 

Property  3.1 

If  x(s)  and  y(s)  are  stable  and  coprime,  then  there  exist  stable  u(s) 
and  v(s)  such  that 

u(s)x(s)  +  v(s)y(s)  =  1  <3.2> 

Also,  if  <3.2>  is  valid  for  stable  x(s)  and  y(s),  then  they  are  coprime. 


Property  3.2 


Xjs) 


<3. 3> 


where  x^(s)  and  y^(s)  are  both  stable,  but  may  not  be  coprime.  Also,  let 


n^(s) 


<3.4> 


be  a  coprime  stable  fractional  representation  for  the  same  r(s).  Then 
there  exists  a  stable  k(s)  such  that 
x^(s)  =  n^(s)k(s) 


y^is)  =  d^{s)k(s) 


<3.5> 


Theorem  3.1  -  Stabilization  Theorom 

For  the  feedback  loop  in  Figure  2  with  transfer  function  ^Cs). 

and  therefore  also  for  the  equivalent  feedback  loop  in  Figure  4  with 

transfer  function  H  (s),  let  the  plant  have  a  coprime  fractional 

''2n‘^ln 

representation 

n  (s) 

P'®’  "  TOT 


There  then  exists  stable  Up(s)  and  Vp(s)  such  that 

Up(s)np(s)  +  Vp(s)dp(s)  =  1 


<3.7> 


Then  for  any  stable  w(s)  such  that  w(s)np(s)  +  v^Cs)  is  not  identically 

zero,  the  compensator 

{-w(s)d„(s)  +  u_(s)}  n_(s) 


<3.9> 


representation  on 

{n  {s)n  (s)} 

Conversely,  every  such  stabilizing  compensator  is  of  this  form  for 
some  stable  w(s).  Note  that  Equation  3.9  guarantees  that  there  are 
no  pole-zero  cancellations  between  the  plant  and  the  compensator. 

Property  3.3 

For  the  same  feedback  loop,  if  c(s)  is  described  by  Equation  3.8 

then 

u  ^  -w(s)n  (s)d  (s)  +  u  (s)n  (s)  <3.10> 

Vl  ''2n“ln  P  P  P  P 

3.3  Problem  Formulation 

Both  the  feedback  loop  and  the  entire  system  in  Figure  4  must 
be  stabilized.  Theorom  3.1  describes  the  parameterization  stabilizing 
the  feedback  loop  portion,  so  the  desired  compensator  stabilizing  both 
quantities  must  be  a  subset  of  this  form.  Referring  to  Figure  4,  define 

n 

and  invoking  Equation  3.10 

^n^^^  ^  ^n^^^  {-Wn(s)np(s)dp(s)  +  Up(s)np(s)}  -  r^(s)  <3.12> 

where  t^(s)  and  '"^(s)  are  described  per  Equation  2.30.  Thus,  the 
objective  here  is  to  stabilize  2^is)  whenever  w^(s)  is  stable.  A 
parameterization  that  stabilizes  Zf,(s)  guarantees  that  v^CO^)  =  y^. 

It  may  not  necessarily  guarantee  that  any  other  derivative  constraints 
are  met,  however.  This  problem  will  be  considered  later  in  this  chapter. 
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3.4  System  Properties  and  Definitions 
In  this  section  groundwork  is  laid  in  preparation  for  analysis 
of  the  system  in  Figure  4  and  Equation  3.12,  leading  up  to  the  IRT 
and  IRP.  Several  properties  and  definitions  will  be  required. 

From  Equation  3.12,  define 


nt(s) 


<3.13> 


=  3757 

n 

to  be  stable  coprime  fractional  representations  of  t^(s)  and  >^p(s). 

As  such,  by  Property  3.1  there  exist  stable  u's  and  v's  such  that 

u^(s)n^(s)  +  v^(s)d.(s)  =  1  <3.14> 

^n  ''n  '^n 

and 

u  (s)n  (s)  +  V  (s)d  (s)  =  1 
n  n  n  n 

Then,  based  on  the  assumption  that  the  plant  is  a  proper  rational 
function,  the  following  property  holds. 


Property  3.4 

The  plant  and  r^(s)  are  such  that 

(d.(s),d  fs))  =  1 
n  ^ 

and  therefore  there  exist  stable  u  (s)  and  v  (s)  such  that 

U  (S)d Js)  t  V  (s)d  (S)  .  1 

n  n  n  ^ 


<3. 15> 


<3.16> 


The  plant  is  assumed  to  be  proper,  so  that  r{p(s)}  ^0,  Thus, 
the  hurwitz  polynomial  which  would  be  used  as  the  common  divisor  of 
both  the  numerator  and  denominator  to  create  np(s)  and  dp(s)  must  be 
of  the  same  degree  as  the  denominator.  Therefore,  r{dp(s)}  =  0.  Such 
a  function  could  have  zeros  only  at  finite  values  of  s. 


On  the  other  hand,  was  described  in  Equation  2.39  to  be  a 

polynomial  in  s,  such  that,  considering  that  some  y^. 's  could  be  equal  to 
zero 

-n  _<  r{r^(s)}  £  0  <3.17> 

If  r{r^{s)}  =  0,  then  r^(s)  is  identically  zero,  so  n^(s)  must  equal 

n 


zero.  Therefore, 


V(s)d  (s)  =  1 
n  n 


<3.18> 


in  which  case  d  (s)  is  proven  coprime  with  d  (s)  by  letting  u  (s) 


equal  zero  and  v  (s)  equal  v  (s).  Also,  if  r{r  (s)}  <  0,  then 


obviously  r{d^(s)}>  0,  and  d^(s)  has  a  constant  as  it  numerator.  In  this 
n  n 


case,  d  (s)  must  be  zero  only  at  infinite  values  of  s.  Thus,  d  (s)  and 
"n  "n 


dp(s)  never  have  common  RHP  zeros,  and  therefore  are  coprime. 
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Definition  3.4 


Define  such  that 


n  (s) 

a,(s)  =  ^ 


n 


n  (s) 
^n 

djil 


u  (s)n  (s)  +  V  (s)d  (a) 
n  n  n  n 


=  1 


<3.19> 


This  function  has  the  same  finite  zeros  as  the  plant,  and  the  same 


finite  poles  as  t^(s),  less  those  that  cancel. 


Definition  3.5 


Define  b^{s)  such  that 


bn(s) 


dJs) 

'n 


dh(s) 


Ub(s)n.  (s)  +  V.  (s)d.  (s)  =  1 
n  n  n  n 


<3. 


This  function  has  the  poles  of  p(s)  as  its  finite  zeros,  and  the  same 
finite  poles  as  t^(s),  less  those  that  cancel. 


The  above  definitions  will  be  needed  later  in  this  chapter,  and 
lead  to  the  following  two  properties. 


Property  3.5 

Invoking  Property  3.2  into  Equations  3.19  and  3.20  leads  to  the 


conclusion  that  there  exist  stable  m  (s)  and  m.  (s)  which  satisfy  the 

^n  °n 

following. 


ina(s)n  (s)  =  n  (s) 
n  n  ^ 


m  (s)d  (s)  =  d.(5) 
n  n  n 


mb(s)nb(s)  =  d  (s)  ;  mb(s)db(s)  =  d^(s) 


<3. 


Property  3.6 

The  functions  mg^(s)  and  mb(s)  are  coprime, 
n  n 

Proof 

Per  Equation  3.7 

Up(s)np(s)  +  Vp(s)dp(s)  =  1 

Then,  substituting  for  n  (s)  and  d  (s)  from  Equation  3.21  yields 
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u  (s)m  (s)n  (s)  +  v  (s)m.  (s)n.  (s)  =  1 

P  P  \  Pn 

{u  (s)n  (s)}m  (s)  +  {v  (s)n.  (s)}m.  (s)  =  1 
^  n  n  ^  “n  n 

which  proves  ni  (s)  and  m.  (s)  are  coprime  by  Property  3.1. 


<3.22> 


Definition  3.6 


Define  d^(s)  as  follows. 


%(s) 

"  FliJ 


<3.23> 


The  function  d^(s)  is  stable. 


Proof 


By  definition 
%(s) 
”  mis) 


Multiplying  by  one,  in  the  form  of  Equation  3.22  leads  to 

Vs) 


n  n 


=  dj,(s)u  (s)n3(s)  + 


db(s)v  (s)nb(s)mb(s) 


n  n 


<3.24> 


niu(s)d.(s)  d.(s) 
n  n  _  n 
m  (s)d  (s)  "  d  (s) 
n  n 


%(s)  dg(s) 
n  _ 

m^JsT  “  djj(s) 
n 

Then,  substituting  Equation  3.25  into  Equation  3.24  yield 

db(s)vD(s)n^(s)da(s) 


<3.25> 


dj,(s)  =  djj(s)u  (s)n^(s)  +  — - ^ 


n  n 


d  (s)  =  d.  (s)u  (s)n  (s)  +  v  {s)n.  (s)d  (s) 
n  n  ^  n  ^  n  n 

which  is  stable  since  it  is  the  product  and  sum  of  stable  functions. 


Note  that  qualitatively  d  (s)  can  be  thought  of  having  as  its 

finite  zeros  the  finite  poles  of  t^(s)  that  are  not  common  with  the 
finite  poles  and  zeros  of  the  plant.  Two  more  properties  will  be 
needed  concerning  d  (s)  before  the  preceeding  definitions  and  properties 


are  used  to  analyze  Equation  3.12. 


Property  3.8 

d*(s)  =  m,{s)m.  {s)d.{s) 


<3.26 


5E 


From  Equation  3.21 


m.  (s)d  (s)  =  d.(s) 
n  °n  n 


<3.27> 


and  from  Definition  3.6 


or 


d^(s)m^(s)  =  djj(s) 


<3.28> 


’n  ~n  n 

Therefore,  substituting  d.  (s)  in  Equation  3.28  into  Equation  3.27  yields 


m^(s)mjj(s)d^(s)  =  d^(s) 


EOP 


Property  3.9 

The  function  d  (s)  is  coprime  with  the  quantity  ng(s)n|j(s)n^(s) , 


n 


so  that  there  exists  stable  u  (s)  and  v  (s)  where 

^n  n 


Ug(s)ng(s)njj(s)n^(s)  +  v^(s)d^(s)  =  1 


<3.29: 


Proof 

This  proof  is  in  four  stages.  First,  (d  (s),n  (s))  =  1  is  proved 

n  ^n 

as  follows. 


u,.(s)n.(s)  +  v.(s)d^.(s)  =  1 


and  invoking  Property  3.8  leads  to 


{u^(s)}n.(s)  +  {v.(s)m  (s)m.  (s)}d  (s)  =  1 
n  n  n  n  n  n 

Next,  (d  (s),  n. (s))  =  1  is  shown  as  follows. 

‘'n  °n 

Uh(s)n.  (s)  +  V.  (s)d.  (s)  =  1 
n  n  n  n 

and  substituting  for  d.  (s)  from  Equation  3.21  yields 

°n 

dt(s) 

U|,(s)nJs)  +  vjs)  =  1 

n  n  “n 

Then,  substituting  for  d.(s)  according  to  Equation  3.26  leads  to 

^n 

{iiu(s)}nb(s)  +  {v.(s)m  (s)}d  (s)  =  1 

n  n  n  ‘^n  n 

Thirdly,  {d_(s),  n,(s))  =  1  can  be  shown  as  follows. 

^n  ®n 

li  (s)n  (s)  +  V  (s)d  (s)  =  1 

n  n  n  n 

and  substituting  for  d^(s)  according  to  Equation  3.21  and  invoking 

®n 

Equation  3.26  leads  to 

dt(s) 

n 

{u  (s)}n  (s)  +  {V  (s)m.  (s)}d  (s)  =  1 

n  n  n  n  n 

Finally,  multiplying  the  three  preceeding  results  together  leads  to 
the  following  result.  Note  that  to  save  space,  the  variable  s  and 
subscript  n  have  been  omitted  from  the  equations. 


k„{s)  = 


Up(s)np(s)nt(s)  -  x^{s) 


<3.34 


are  stable,  where 


n„(s)d.(s) 

If  ^-1  -  "  " 

— urr 

n 


and  where  x„(s)  is  stable, 
n 


<3.3£ 


Proof 


Equation  3.12  is  rearranged  as  follows. 

z„(s)  =  t^(s){-w^(s)np(s)dp(s)  +  Up(s)np(s)}  -  r^^Cs) 


n^(s) 


“  3757  ■  3117 

^n  '"n 

z^(s)d^(s)d^(s)  =  n^(s)d^  (s){-w^(s)n  (s)d  (s)  +  u  (s)n  (s)} 


n  n 


n  n 


n  n 


z^(s)d^(s)d^(s)  +  w^(s)np(s)dp(s)n^  (s)d^(s)  -  Up(s)np(s)n^(s)d^(s) 


=  -r>p(s)d^(s) 
n  n  ^ 


n  n 


n  n 


zjs)cljs)  +  w„(s)n„(s)d„{s)n.(s)  -  u„(s)n„{s)n^{s)  = 


■n^{s)d^(s) 
n  ^n 


x^(s)  must  be  stable  since  it  is  equal  to  the  product  and  sum  of  stable 
functions  when  z^(s)  and  w^(s)  are  stable.  Then,  substituting  Equations 
3.26  and  3,21  into  the  above  equation  and  rearranging  leads  to  (omitting 
the  variables  and  subscript  n) 

zm,m.  d^  +  wm,n,m.  n.  n.  -  u„n„n.  =  -x 
abc  aabbt  ppt 

(zd^  +  wn,n.  n.)m,m.  =  u„n„n.  -  x 
c  aot  ao  ppt 


u_n_n. 


zd^  +  wn,n.  n.  = 
c  a  b  t 


<3.36 


If  Zp(s)  and  w^(s)  are  stable,  then  the  left  side,  and  therefore  the 
right  side  of  Equation  3.36  is  stable.  Also,  in  this  case,  if  Equation 


3.36  is  multiplied  by  either  m  (s)  or  m.  (s) 

n  °n 

u  n  n.  -  X 

(zd^  +  wn,n.  n.  )m,  =  -£..£..1 - -  k 

'  c  a  b  t  a  m. 


<3.37 


u  n  n.  -  X 

(zd  +  wn,n.  n.)m.  =  -P- - =  g  <3.3£ 

c  a  b  t  b  m  ^ 

Q 

the  left  side  of  the  equation  is  still  stable.  Therefore  kp(s)  and 
g^(s)  are  stable.  Thus  far  we  have  proven  that  in  order  for  a  solution 
to  exist,  the  right  side  of  Equations  3.36,  3.37,  and  3.38  must  be 
stable.  To  prove  that  stable  kp(s)  and  g^Cs)  is  sufficient  for  the 
existence  of  stable  2p(s)  and  w(s),  first  multiply  the  right  side  of 
Equation  3.36  by  one,  in  the  form  of  Equation  3.29. 


zd^  +  wn,n.  n. 
c  a  b  t 


-  X 

ppt 


(u  n  n  .  x)v  -  “>“3 

={  -P-  -P- ■  ■ - ^}dc  +  B-g-j - ^}n  n,n. 


Therefore,  if  Equation  3.36  is  stable,  then  particular  stable  solutions 


obviously  exist  for  2^(s)  and  w^(s).  Also,  <3.36>  is  stable  whenever 
g^(s)  and  kp(s)  are  stable.  This  can  be  proven  by  multiplying  by 
Equation  3.22,  as  follows. 


u_n_n^  -  X  u_n_n^  -  x 


{u  n^m  +  v^n.m.  } 
p  a  a  p  b  b 


u_n_n^  -  X 


u_n  n^  -  X 


P  3 


=  ku  n,  +  gv  n. 
pa  p  b 


Thus,  if  k^(s)  and  g^(s)  are  stable,  the  right  side  of  Equation  3.29 
is  stable,  and  therefore,  also  the  left.  Finally,  to  show  that  k^(s) 
and  gj,(s)  must  be  stable  if  3.36  is  stable,  let 


u  (s)n  (s)n^(s)  -  x^(s) 


m  UJm.  (sj 
n  n 


where  z’(s)  stable.  Then  multiplying  by  m^(s)  yields 
"  ®n 

k  (s)  =  m  {s)z'(s)  <3. 

n  a^  n 

which  must  be  stable  if  Zp(s)  is  stable.  Also,  multiplying  <3.40>  by 
mjj(s)  leads  to 


9n<  =  >  ’ 

n 


which  is  stable.  Therefore,  it  is  both  necessary  and  sufficient  that 
x„{s),  k„(s)  and  g„(s)  are  stable  in  order  for  z„(s)  and  w„(s)  to  both 


be  stable. 
EOP 


The  above  property  defines  the  basic  mathematical  commodities  that 
must  be  stable  in  order  for  a  solution  to  Equation  3.12  to  exist. 
However,  the  commodities  are  derived  mathematically  from  the  original 
system,  and  give  a  designer  no  "feel"  as  to  their  meaning.  In  practice 
one  would  like  the  necessary  and  sufficient  conditions  to  be  related 
simply  to  basic  system  properties  which  are  not  derived  mathematically. 
With  this  goal  in  mind,  the  quantities  given  in  Property  3.10  are 
further  analyzed  as  follows. 

Property  3.11 
x„(s) 

mjj(s) 

n 

is  stable. 

Proof 

Again,  omitting  the  variable  s  and  subscript  n  for  compactness 


by  definition,  so  that 
x_  .  "r‘‘t 

Then,  using  Property  3.8  yields 


n  m,m.  d 


n  m.d 


Multiplying  by  one  in  the  form  of  Equation  3.16  leads  to 


rrijj  dj,  rp  r  rp  p 

n  m  d^v^„d„ 

=  n  m  d  u  +  ^  ^  P 

r  a  c  rp  d^ 

Then,  substituting  for  dp(s)  according  to  Equation  3.21  leads  to 

-  "A^cVp  * - — 

n  m,m.  d_ 

=  n  m  d  u  + - T-P-^  V  „n. 

r  a  c  rp  d  rp  b 


—  =  nmdu„  +  xv  „n, 
m^j  r  a  c  rp  rp  b 


^-1 


<3.43> 


The  right  side  of  Equation  3.43  is  stable,  and  therefore  so  is  the  left. 
EOP 


Property  3.11  in  itself  provides  no  further  immediate  insights, 
but  it  does  lead  directly  to  the  following  lemma  which  does. 

Lemma  3.1 

Equation  3.34  admits  stable  solutions  iff 

(d.{s),d  (s))  =  1  <3.44> 

If  this  is  the  case,  then  there  exist  stable  u^  (s)  and  (s) 

nP  nP 


such  that 

u.  _(s)d.(s)  +  V.  n(s)d„(s)  =  1 


<3.45: 


Rearranging  Equation  3.34  as  follows 
Li.,(s)n_(s)n^(s)  -  x_(s) 


u..(s)n^{s)n*  s) 


u„(s)n„(s)n*(s) 


If  l<^(s)  is  to  be  stable,  then  the  left  side  of  <3.47>  is  stable  by 
invoking  Property  3.11.  Therefore,  the  right  side  must  be  stable  also. 
Stability  of  the  right  side  of  <3.47>  sufficiently  guarantees  stability 
of  k^(s)»  since  in  this  case  Equation  3.46  represents  k^(s)  as  the  sum 
of  two  stable  functions. 


Referring  to  Equation  3.44,  if  either  d^(s)  or  d  (s)  are  miniphase, 

then  they  are  obviously  coprime.  This  can  be  proven  by  multiplying  the 
mini  phase  function  by  its  inverse  (yielding  1)  and  the  other  function 
by  zero.  Thus,  we  have  narrowed  the  proof  down  to  the  case  where  both 
functions  have  zeros  in  the  RHP. 


Equation  3.20  effectively  cancels  the  common  zeros  between  dp(s) 

and  d.(s)  in  creating  n.  (s)  and  d.  (s),  which  in  effect  are  then  equal 
^n  °n  °n 

to  the  first  two  functions  less  the  common  zeros.  Then,  by  Equation  3.21, 


s 


Thus,  m.  (s)  is  the  original  function  divided  by  itself  without  the 
°n 


common  zeros,  yielding  a  function  only  containing  all  zeros  that  are 


common  between  dp(s)  and  d^(s). 


If  d.(s)  and  d  (s)  are  not  coprime,  then  the  above  argument  dictates 


that  m.  (s)  has  RHP  zeros,  and  is  not  miniphase.  Then,  to  stabilize 
<3.47>,  it  is  necessary  that  (m.  iu_n  n.  ).  But,  d  (s)  is  coprime  with 


Up(s)np(s),  because 


{l}Up{s)np{s)  +  {Vp{s)}  dp(s)  =  1 


Thus,  mjj(s)  is  also  coprime  since  it  only  has  zeros  from  d  (s), 


Therefore,  m.  (s)  must  divide  n^(s),  and  in  order  for  this  to  occur,  it 
“n  ^n 


must  have  common  zeros.  But,  m.  (s)  also  derived  its  zeros  strictly 


from  d.(s)  which  has  no  common  RHP  zeros  with  n.(s).  Therefore,  m.  (s) 
^n  ’'n  ^n 

can  not  divide  n.(s).  This  proves  that  m.  (s)  must  be  miniphase,  also 
S  °n 

proving  that  d  (s)  and  d^(s)  must  be  coprime. 


EOP 


Lemma  3.1  provides  the  necessary  and  sufficient  condition  in 
order  to  stabilize  Equation  3.34.  It  is  interesting  to  note  that  the 
lemma  simply  requires  that  the  plant  not  have  any  common  RHP  poles 
with  the  input,  less  perhaps  up  to  (n  +  2)  poles  at  s  =  0  of  the  input. 
For  most  systems,  this  of  course  should  present  no  difficulty,  and 
therefore  is  not  a  very  strict  condition.  The  condition  necessary  to 


stabilize  Equations  3.33  and  3.35  must  still  be  derived.  The  following 
property  leads  towards  Lemma  2,  which  does  this. 


Property  3.12 

If  k  (s)  is  stable,  then  g„(s)  is  stable  iff  d^(s)  divides  d  (s), 

"  '^n  ^n 

in  which  case  x^(s)  is  also  stable. 


Proof 


must  then  be  focused  on  the  case  where  it  is  not  mini  phase.  Then, 


using  Equation  3.21 

u„(s)n„(s)n. (s)  u  (s)m  (s)n  (s)n.(s) 

P  P  .  n  _  P  %  S 

inTT?)  m  As) 

n  n 

=  u  (s)n  (s)n  (s) 

^  n  n 

which  is  stable.  Thus,  rearranging  <3.33>, 

Up(s)np(s)n^(s)  -  x^(s) 

^n^^^  ^  m  (s) 

^n 

X  (s) 

=  u  (s)n3(s)n^(s)  - 
^  n  n  a^ 

Xn(s) 

gn(s)  +  Up{s)n^(s)n^(s)  =  -  ^ 

n 

Equation  3.49  shows  that  the  right  side  must  be  stable  if  9^(s)»  and 
therefore  the  left  side,  is  stable.  Equation  3.48  shows  that  if 


<3.' 


<3.- 


is  stable,  then  so  is  g^(s).  Thus,  stabilizing  this  term  is  necessary 
and  sufficient  to  stabilize  g^(s).  Then,  invoking  Equations  3.35  and 
3.26  leads  to 

x„(s) 

n  _  _  _  n  n 

m,(s)  ■  d  {s)m  (s) 
n  n  n 

n  (s)m  (s)m.  (s)d  (s) 

=  n  n  n  n 

cljs)m  (s) 

•^n  ^n 


n  (s)m.  (s)d  (s) 
n  n  n 


<3.50> 


which  must  be  stable.  Analyzing  <3.50>,  one  notes  that  d  (s)  must 
divide  the  numerator,  since  we  are  concerned  with  d  (s)  not  miniphase. 

But,  it  can  not  divide  n  (s)  even  partially  since  they  are  coprime. 

Also,  in  the  case  that  k^(s)  is  stable,  during  the  proof  of  Lemma  3.1, 

it  was  shown  that  m.  (s)  is  miniphase,  so  d  (s)  cannot  divide  m.  (s)  either. 

°n  ^n  “n 

Therefore,  in  order  for  g_(s)  to  be  stable  when  k_(s)  is  stable,  d  (s) 

n  ”  ”n 

must  divide  d^(s).  Note  that  this  condition  also  guarantees  x^(s)  is 
n 


stable,  since 


<l^(s) 


x„(s)  =  tl^(s)m^(s)n.^(s) 


n  n  n  r: 


The  above  property  still  does  not  relate  directly  to  system 


quantities,  since  d  (s)  is  derived.  Thus,  the  following  lemmas  is  given. 


whose  proof  further  analyzes  the  condition  in  Property  3.12  in  order  to 
determine  what  is  physically  meant  by  the  statement. 


Lemma  3.2 

If  stable,  then  9^,(5)  and  x^(s)  are  both  stable  iff 

r{t(s)}  +  r{p(s)}  £  M  +  1  <3.51: 

where  M  is  the  lowest  order  initial  derivative  constraint  not  equal  to 
zero,  i.e.,  y^  =  0  for  i  <  M,  and  yj^  i-  0. 

Proof 

If  l<f,(s)  is  stable,  then  by  Property  3.12  gj^(s)  and  x^(s)  are 

stable  iff  (d  |d  ).  Thus,  this  condition  must  be  analyzed  to  see 
^n  ^n 

when  it  occurs. 


Before  proceeding  further  however,  the  following  shorthand  notation 
is  introduced  which  will  prove  useful.  First,  let  0^(n)  represent  a 
stable  (hurwitz)  polynomial  of  order  n,  and  let  Oy(n)  represent  a 
general  polynomial  of  order  n,  where  n  is  an  integer.  In  this  notation 
if  n  is  less  than  zero,  then 


0(n)  = 


OT^ 


<3.52 


Equation  2.39  defined 


r  Js)  =  E 


i  =  0 


and  t^(s)  =  s'^''‘^t{s) 


Let  M  be  the  lowest  degree  initial  response  constraint  not  equal  to 
zero,  so  that  y^.  is  not  zero,  but  {yf,,...,yf.  ,}  =  0. 


45 


Then, 


r^(s)  =  +  ...  +  =  Oj^{n-M+l) 

The  copritne  representation  of  r^(s)  is  then 


<3.53> 


Oy(n-M+l) 

Ojn-M+1) 

r„(s)  =  Ojn-M+1)  =  - 

Og(n-M+l) 


n^(s) 

n 

n 


<3.54> 


The  plant  is  such  that  r{p(s)}  ^  0  by  assumption.  Therefore,  the 
coprime  representation  of  p(s)  is  developed  as  follows. 


u'-^p' 


<  y„ 

p  -  j'p 


W 

370^  =  3^  •  ’'pi)', 


<3.55> 


For  the  coprime  representation  of  t^(s)  =  s'’  '^t(s),  it  is  necessary 
to  know  how  many  poles  at  s  =  0  t(s)  has.  Thus,  let 


t(s)  = 


0„(xJ 


<3.56> 


where  c  is  a  non-negative  integer,  and  where  Oy(y^)  represents  a 
polynomial  that  is  not  divisable  by  s.  Also,  if  c  is  positive,  then 
Ou(Xt)  is  not  divisable  by  s.  Then, 


t.(s)  = 


W" 


0  (A)0  (B2) 


<3.63 


To  calculate  d  (s).  Equation  3.26  is  used,  so  that 
dt(s) 
n  n 

Earlier  in  this  paper  it  was  shown  that  m.  (s)  must  be  miniphase  in  order 

°n 

for  a  solution  to  exist.  Thus,  let  x.  =  y.  ,  so  that 


K] 


w  ww  w 

w  ■ 


Finally,  d  (s)  must  divide  d  (s),  so  that 

n 


0  (y  -A)0  (y  )0s(Bi)0s(yb)0s{n-M+l) 


0s(Bi)0s{B2)0s{Xb) 

n 

must  be  stable.  Note  that  the  denominator  contains  only  LHP  poles, 
since  it  is  the  product  of  stable  polynomials.  Therefore,  the  quotient 
is  stable  iff  the  degree  of  the  denominator  is  equal  to  or  greater  than 
the  degree  of  the  numerator.  This  requires  that 
Bi+B2+Xb  >.  n-M+l+y^-A+yp+Bj+yb 
or,  cancelling  like  terms  and  noting  that  yb  =  Xb, 

B2  >  n-M-l+2+Y^-A+yp 


M+1  >  yb‘''yp“A-B2+n+2  =  B^ 

Then  analyzing  <3.67>  further, 

=  y^+yp-A+n+2-max{Xp-A+max{y^,n+2-c+x^},y^-A+yp} 

=  (yb+yp-A+n+2)+min{-Xp+A-max{y^,n+2-c+x^},-y^+A-yp} 
since  -max{a,b}  =  min{-a,-b}.  Then,  adding  the  term  inside  the 


parenthesis  to  both  sides  inside  the  min  function, 

Bj  =  min{y^+yp-A+n+2-Xp+A-max{y^,n+2-c+x^},  y^+yp-A+n+2-y^+A-yp} 

=  mi  n{y^+yp+n+2-Xp-max{y^ ,n+2-c+x^} , n+2} 

=  min{(y^+yp+n+2-Xp)+min{-y^,-n-2+c-x^},n+2} 

=  min{min{y^+yp+n+2-Xp-y^,y^+yp+n+2-Xp-n-2+c-x^},n+2} 

=  min{min{y  -X  +n+2,y.+c-x.+y  -X  },n+2}  <3 

P  P  ^  ^  P  P 

The  first  quantity  in  the  interior  min  function  is  yp“Xp+n+2.  But, 
according  to  Equation  3.55,  Xp_<yp.  Therefore, 

y„-x.,+n+2  >  n+2 
P  P  - 

so  <3.68>  can  be  further  reduced  to 


<3.68> 


B^  =  min{n+2,y^+c-Xj+yp-Xp} 

M+1  i  min{n+2,y^+c-x^+yp-Xp} 


<3.69> 


Thus,  in  order  for  d  (s)  to  divide  d  (s),  the  inequality  in  Equation  3.69 

n  n 

must  be  met.  But,  M  _<  n,  so  M+1  cannot  be  equal  to  or  greater  than  n+2. 
Therefore,  the  inequality  must  relate  to  the  second  term  within  the  min 
function.  Also,  noting  that 
yt+c-^t  ^  '"{t(s)} 


yp-Xp  -  r{p(s)} 

we  can  finally  state  that  the  necessary  and  sufficient  condition  is  that 


m+1  >  r{t(s)}  +  r{p(s)} 


<3.70> 


The  above  inequality  is  based  on  dj.(s)  according  to  Equation  3.57, 


leaving  the  inequality  based  on  d.{s)  according  to  <3,58>  yet  to  be 

proved.  This  proof  is  similar  to  the  one  just  given,  and  is  also 
lengthy,  so  in  order  to  condense  it  many  of  the  comments  between  steps 
will  be  omitted.  First,  however,  it  will  prove  useful  to  modify  <3.58> 
so  that 


d.(s) 

^n 


0^(max{x^,c-n-2+y^}j 


;  c-n-2  >  0 


Then, 


n  (s)  )  0  (max{x  ,c-n-2+y  }) 


Ou(Xp-A' )0^(max{Xt,c-n-2+yt}) 
0^^(c-n-2+y^-A-)0^(yp) 


n  (s) 
^n 


0,(x  -A')0s(B4) 


W 


where 


84  =  max{x^,c-n-2+y^} 

Bg  =  max{Xp-A'+B4,c-n-2+y^-A'+yp} 


d.(s) 

^n 


0y(c-n-2+y^) 

■“W” 


n„{s)  0.,{x„) 


OjBc) 


0..{A')0  fBj 


for  a  solution  to  exist.  Then, 

=  c+y^+y  -Xp'*fflin{-x^,-c+n+2-y^} 

=  min{c+yp+yp-Xp-Xj.c+yp+yp-Xp-c+n+2-y^) 

=  iiHn{yp+c-Xptyp.Xp,yp-Xptnt2) 

Also,  y  >  X  ,  and  n  ^  M,  so  M+1  can  not  be  greater  than  yn"AD'*'h'*'2- 

nr'  r  n 

Thus,  stability  requires  that 

H+l  >  (y^tc-Xp)  t  (yp-Xp) 

Finally, 

M+1  ^  r{t(s)}  +  r{p(s)}  <3.81> 

Equations  3.70  and  3.81  are  equivalent  to  Equation  3.51,  thus  proving 
Lemma  3.2. 
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If  the  conditions  are  met  in  Lemma  3.1  and  Lemma  3.2,  then  a 
compensator  exists  that  will  stabilize  the  feedback  loop  and  meet  the 
initial  n-th  derivative  condition.  At  this  point,  however,  there  is 
no  guarantee  that  such  a  compensator  will  meet  all  other  constraints 
also,  which  is  the  final  goal.  Thus,  a  complete  parameterization  meeting 
the  n-th  derivative  constraint  will  now  be  derived.  Then  this 
parameterization  will  be  analyzed  to  see  when  it  will  also  meet  the 
constraints  tYg* • • • ’^n-l^' 

Lemma  3.3 

For  the  system  in  Figure  4,  described  by  Equation  3.12,  iff  the 


conditions  r{t(s)}  +  r{p(s)}  £  M+1,  and  (d.(s),d  (s))  =  1  are  met, 

then  a  compensator  exists  that  stabilizes  both  w(s)  and  Zp(s).  In 
this  case,  c(s)  is  described  by  Equation  3.8,  where  the  complete 
parameterization  for  w(s)  is 


da(s) 


w(s)  =  w^(s)  =  (Up(s)na(s)n^(s)-n^{s)  j-^Yy^u^(s)  +  da(s)e(s)  <3.82> 


whenever  w(s)np(s)  +  Vp(s)  is  not  equal  to  zero.  The  function  e(s)  is 


arbitrarily  stable,  a„(s)  is  per  Equation  3.19,  and  u.(s)  (and  v.(s)) 


is  stable,  and  satisfies 


U5(s)n  (s)d  (s)n^(s)  +  v^(s)da(s)  =  1 


<3.83> 


With  w^(s)  described  per  Equation  3.82 


d  (s) 


z„(s)  =  (Up(s)nj(s)nj(s)-n^(s)  v^ls) 


n  ^  n 


<3.84> 


Proof 

To  admit  stable  solutions.  Property  3.10  requires  stable  g^(s), 

k  (s),  and  x  (s).  By  Lemma  3.1,  (d.(s),d  (s))  =  1  is  necessary  and 
n  n  ^n  ^ 

sufficient  to  stabilize  k^(s).  By  Lemma  3.2,  r{t(s)}  +  r{p(s)}  j<  M+1 
is  N  &  S  to  stabilize  g_(s)  and  x  (s),  by  guaranteeing  that  (d  |d  ). 

n  fi  I  _ 

n  n 

Thus,  the  conditions  are  justified. 


Assume  these  conditions  are  met.  Then,  Equation  3.12  has  been 
modified  to  Equation  3.36  (omitting  the  variable  s  and  subscript  n  during 


manipulation) , 


u_n_n.-x 


zd^  +  wn,n.  n.  = 
C  d  D  L 


Vo"t- 


Also,  b^(s)  was  defined  by 


d.(s)  "b<=) 


But,  (d  (s),d.(s))  =  1,  so  we  can  let 
P  ^n 


d  (s)  =  n.  (s)  ;  m.  (s)  =  1  ;  d^(s) 

^  n  n  n 


d,(s) 


Then,  by  Equation  3.26, 

d.(s)  =  m  (s)m.  (s)d  (s) 
n  n  n  n 


'  m  (s)d  (s) 
n  n 


But,  by  <3.21>, 

d.(s)  =  '"a(s)d  (s) 
n  n  n 


so  that 

m  (s)d  (s)  =  m  (s)d  (s) 
n  n  n  n 


d  (s)  »  d  (5) 

n  n 


Then,  by  Property  3.9,  there  exist  stable  u,{s)  and  v„(s)  such  that 


U  (s)n  (s)n.  (s)n.(s)  +  v  (s)d  (s)  =  1 
n  n  n  b  n  n 

This  can  be  modified  to 

^5(5)0  {s)d  (s)n.(s)  +  V  (s)d  (s)  =  1 
n  n  ^  n  n  n 

Also,  <3.85>  can  be  modified  as  follows. 


Li_n_n.  -  d„ 


zd^  +  wn,n.  n.  = 
c  a  b  t 


zd,  +  wn,d„n.  = 
a  apt 


“pVt-  “r 


_  P  P  t 

"r'^t 

m. 

a 

a  r 

-  P  P  t 

"r"a'‘a 

n 

m  d 

_E 

a  r 

"a 

zd,  +  wp,d„n,  =  u„n,n,  -  n_  <3.8 

a  apt  pat  rd^ 

The  right  side  is  stable  since  dr^ld^  (d^  =  d^  ). 

Thus  to  find  particular  solutions,  multiply  by  Equation  3.87  as  follows 


zd^  +  wn^dn.  =  (u  n^n.  -  n^  ■j^)-(v^d^  +  u^n^d  n.) 
a  apt  pat  r  d^'  '  s  a  s  a  p  t 


d,  d 

'  ‘(“p"a"t  -  "r  *  ''“pVt  -  "r  3:>“s>"aVt 


Particular  solutions  are  then  equal  to 


d,(s) 


wj(s)  =  (Up(s)rtj(s)nj(s)  -  n^Js) 


zj(s)  =  (u  (s)nj(s)n^(s)  -  n^(s)  g-fjy)»5(s) 


<3.90 


n  n 


To  find  the  complete  solution  set,  we  also  need  the  homogeneous 
solutions.  Typcially,  this  is  accomplished  by  guessing  the  homogeneous 
solutions  that  cause  the  left  side  of  the  equation  (here.  Equation  3.89) 
to  equal  zero.  Then,  the  guesses  are  tested  to  verify  that  they  are  in 
fact  the  homogeneous  solutions.  Thus,  first  guess  that 


w"  =  d  e 

a 


z  =  -n^d  n.e 
apt 


<3.91 


where  e(s)  is  arbitrarily  stable,  so  that 

h  h 

z  d^  +  w  n,d  n.  =  (-n,d  n.e)d,  +  d,en,d  n.  =  0 
a  apt  apt  a  aapt 

as  desired,  verifying  that  <3.91>  are  homogeneous  solutions.  We  must 

also  test  that  all  homogeneous  solutions  are  of  the  form  of  <3.91>.  To 

h  h 

do  this,  assume  that  ^(s)  and  zj^(s)  are  stable,  and  satisfy 


z*^d.  +  w^n.d  n.  =  0 
—  a  —apt 

and  define  the  function  e(s)  by 


e  =  ^ 
-  d. 


<3.92 


<3.93 


Clearly, 


57 


and  it  follows  from  <3.89>  that 


.  -ed^n^d  n. 

_h  — aapt_ 

z  - - 3 —  =  -n^d^n.e 

—  d  a  p  t— 


<3.95> 


h  h 

showing  that  w"  and  have  the  form  of  Equation  3.91.  Finally,  to 

insure  the  completeness  of  the  solution  set,  ie(s)  must  be  shown  to  be 

stable.  To  do  this,  the  coprimeness  properties  of  n,(s)  and  d,(s),  and 

n  n 

of  n.(s)  and  d.(s)  are  needed,  as  follows. 

^n  ^n 


£  “  t^Va  *  Va> 

a 


.h..  ,  -HI  a  a/ 


.  w"v,  *  ^Vt  * 


h  w  u^n^v.d. 

_  +  dL..l_.a_t_t  ^  r  .A^Lt  t 


Invoking  m,d,  =  d.  and  (d.  ,dp)  =  1, 
a  a  c 


.  w  u,n,v.m^d,  w.  u,n,a.n. 

1  =  w'»,  *  ~  ^  r  ^  ^  ^ 


h  h  w. u^n^u.n. v.„d^ 

e  =  w^v,  +  w^u.n.v.m,  +  w.  u  n  u.n.u.  m  +  ^  ^  ^  ■■  ■  ■  ■  P 
—  —  a  —  aata  -Hiaattt_a  a 

P  a 


<3.96> 


Also,  by  modifying  Equation  3.92  yields 


u 

■h  ■  ''a‘‘p''t 

- 


<3.97> 


Thus,  substituting  <3.97>  into  <3.96>  leads  to 


h  h  h 

e  =  w"v,  +  w"u,n  v.m,  +  w.  n, u.n.u.  m,  -  z  u,u.v.„ 
—  —  a  —  aata  — hatttpa  —  attp 


<3.98> 


which  proves  e(s)  is  stable  since  it  is  the  product  and  sum  of  stable 
functions.  Therefore,  since  the  complete  solution  set  is  the  sum  of  the 
particular  solution  and  the  homogeneous  solution.  Equation  3.82  describes 
the  general  solution  set  for  w(s),  and  Zp(s)  is  described  by  <3.84>. 

To  complete  the  proof,  we  must  show  that  at  least  one  stable  e(s) 
leads  to  a  w(s)  such  that  w(s)np(s)  +  Vp(s)  is  not  identically  equal  to 
zero,  in  which  case  since  this  is  the  denominator  of  the  compensator, 
the  resulting  compensator  function  would  be  undefined.  Thus,  again 
omitting  the  variable  s  and  subscript  n. 


wn„  +  V  =  {(u  n  n.  -  n  -j^)u  +  d  e}n  +v 
P  p  pat  rd's  a  p  p 


=  ((u  n  n.  -  n„ 

pat  r  s  p  p  a  P 


<3.99> 


The  function  d  (s)  is  not  zero  since  it  is  the  denominator  of  a  (s). 

®n  " 

Also,  if  np(s)  is  not  zero,  then  w{s)np{s)  +  Vp(s)  is  a  non-trivial 
function  of  e(s),  and  is  therefore  not  identically  zero  for  all  e(s). 
If  HpCs)  is  equal  to  zero,  then  <3.7>  implies  that  Vp(s)  is  miniphase, 
and  therefore  can  not  be  equal  to  zero.  Thus,  in  this  case 
w(s)np(s)  +  Vp(s)  =  Vp(s)  is  not  equal  to  zero.  Therefore,  the  proof 
is  complete. 


3.6  Simultaneous  Derivative  Constraints 


Clearly,  Lemma  3.3  could  be  used  to  find  a  compensator  meeting  all 


every  separate  constraint,  from  which  the  resulting  parameterization 
would  be  the  intersecting  parameterization  set  of  all  the  individual 
parameter! zations.  This  could  require  a  great  deal  of  computation 


since  t^(s),r^(s),  and  the  resulting  computed  functions  such  as  ^^(s) 
depend  on  n.  This  normally  will  not  be  required  though.  Usually  (and 
possibly  always)  solving  just  the  n-th  derivative  case  also  will  guarantee 
solving  the  0-th  through  (n-l)th  cases.  A  set  of  necessary  and  sufficient 
conditions  for  this  guarantee  have  not  been  found,  but  the  following 
lemma  proves  that  usually  only  the  n-th  derivative  parameterization  will 
need  to  be  computed. 


Lemma  3.4 

Suppose  the  conditions  are  met  to  apply  Lemma  3.3  for  the  n-th 
derivative  initial  response  constraint,  so  that  w^(s)  is  computed  per 
Equation  3.82,  resulting  in  per  Equation  3.84.  Then  any  compensator 

based  on  Equations  3.82  and  3.8  will  stabilize  not  only  z^(s),  but  also 
Zq(s)  through  z^_j(s),  whenever 

1.  The  input  t{s)  has  fewer  than  three  poles  at  s  =  0.  <3.100> 

or 

2.  If  the  above  condition  is  not  met,  then  <3.101> 

{{#poles  t(s)  at  s  =  0)  -  3}  <  {#  zeros  p(s)  at  s  =  0}. 


Proof 


Zp(s)  is  given  by 

z^(s)  =  t^(s){-w(s)np(s)dp(s)  +  Up(s)np(s)}  -  r^{s) 


=  s""2t(s)H,  ,  (s)  -  E 

''2  1  i=0  ' 


<3.102> 

Where  ,,  (s)  is  stable  whenever  w(s)  is  stable,  as  is  the  case  in 
''2^1 

Equation  3.82.  Then,  manipulating  <3.102>  leads  to 


V-=  s-^t(s)H^  (s)  -  E  s"  -y. 

2  1  i=0 


z„(s) 


-V-  *  =  )  -  f.,  ="”yi 

But,  using  the  same  set  of  constraints. 


z„  .(s)  =  „{s)  - 

n-i  V2UJ  1 


=  s"'^^t{s)H  (s)  -  E 

2  1  i=0 

Comparing  <3.103>  and  <3.104>  leads  directly  to 
z  (s) 

Vl<=)  ’  J'n 


<3.103> 


<3.104> 


<3.105> 


Thus,  since  y^  is  a  constant,  its  sufficient  but  possibly  not  necessary 
to  prove  the  z^_j(s)  is  stable  by  proving  s  divides  z^.  Then,  if  w(s) 
is  given  by  <3.82>,  z^(s)  is  given  by  <3.84>  to  be 

djCs) 


z„(s)  =  {u  (s)n  (s)n  (s)  -  n  (s)  j-^)  v  (s) 

y  “n  '-n  'n  r'  '  " 


n  n 


ng(s)d  (s)  n^{s)e(s) 


where  dr^(s)  divides  d^(s),  so  that  z^{s)  is  stable.  The  above 

n 

equation  can  be  modified  to 


z  (s)  =  n.(s)n  (s){u  (s)v  (s)  -  d  (s)e(s)} 
n  n  n  ^ 

da(s) 

n 


<3.106> 


Clearly,  n^(s)  will  always  be  divisable  by  s,  because  <2.39>  describes 
n 

r^(s)  as  a  polynomial  in  s  with  no  constant  term.  Thus,  z^(s)  will  be 

divisable  by  s  whenever  either  n.(s)  or  n,(s)  are  divisable  by  s 

n  n 

(possibly  other  times  too),  since  z^(s)  will  then  be  the  sum  of  terms 
divisable  by  s. 


Therefore,  first  consider  n.(s),  as  follows. 

^n 


Generally, 


t(s)  = 


s\(y() 


<3.107> 


is  a  representation  for  t(s),  where  the  s  term  represents  the  only 
poles  or  zeros  at  s  =0,  so  that  c  is  any  integer.  This  leads  to 


Obviously,  n^(s)  will  be  divisable  by  s  iff 


<3.108> 


n+2-c  >  0 


<3.109> 


Now,  if  <3.109>  holds  for  n,  then  z^_j(s)  is  stable.  Then,  this  entire 
argument  can  be  repeated  for  the  case  is  stable  to  show  when 

^n-2^^^  is  stable,  which  will  result  in  the  inequality 

(n-l)+2-c  >  0  <3.110> 


Thus,  to  stabilize  Zq(s)  through  2^_2(s),  the  worst  case  inequality, 
which  proves  all  z^.(s),  i  _<  0  £  n-1,  are  stable,  is  the  case  when  n  =  1 
and  (n-1)  =  0.  The  resulting  inequality  is 
1+2-c  >  0 
or 

c  <  3  <3.111> 

which  is  equivalent  to  saying  t(s)  has  fewer  than  three  poles  at  s  =  0. 
Therefore,  <3.100>  is  proved. 


If  n^(s)  is  not  divisable  by  s,  another  condition  that  z  (s)  is 

L  1 1 


divisable  by  s  is  that  n  (s)  is  divisable  by  s.  Pursuing  this  case,  first 


assume  that  n.(s)  is  not  divisable  by  s,  so  that  Equation  3.111  does 
^n 


not  hold,  but  that 


0„(x.) 

tn(s)  =  - 

s"‘"-\{y^) 


<3.112> 


where  c-n-2  ^0  is  a  general  expression  for  t^(s)  based  on  Equation 


3.107. 

Thus, 


d.(s)  = 


0j(max{x^,c-n-2+y^}) 


<3.113> 


Also,  describe  p(s)  generally  by 


5*^0  (xj 
°  0  (y  ? 


<3.114> 

where  D  >  0,  and  0,,(y„)  may  be  divisable  by  s,  so  that  D  represents  the 


number  of  zeros  the  plant  has  at  s  '  0.  Of  course  this  leads  to 


<3.115> 


since  r{p(s)}  ^  0  is  a  general  assumption.  This  of  course  leads  to, 
by  definition. 


5*^0  (x  )  0  (max{x.,c-n-2+y  }) 


<3.116> 


n 


Obviously,  n  (s)  has  max{0,D-c+n+2}  zeros  at  s  equal  to  0.  Thus,  if 
^n 


D-c+n+2  >  0 


<3.117> 


for  all  n  >  1,  then  n,(s)  will  be  divi sable  by  s.  Studying  the  worst 


case  again  of  n  =  1  leads  to  the  general  requirement  that 
D-c+1+2  >  0 


or 

c-3  <  0  <3.118> 

By  our  definitions.  Equation  3.118  is  equivalent  to  the  statement  of 
Equation  3.101,  thus  the  proof  is  complete. 
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3.7  Initial  Response  Theorem  and  Parameterization 
The  exciting  climax  of  this  thesis  is  finally  upon  us!  The 
groundwork  has  all  been  laid,  so  the  Initial  Response  Theorem  and  the 
Initial  Response  Parameterization  can  now  be  given. 

Theorem  3.2  -  Initial  Response  Theorem  (IRT) 

For  the  system  in  Figure  2,  assume  that  r{p(s)}  ^0.  Also  assume 


64 


that 


or 


1.  The  input  t(s)  has  fewer  than  three  poles  at  s  =  0, 


2.  If  the  above  condition  is  not  met,  then 

{(#  poles  t(s)  at  s  =  0)-3}  <  {#  zeros  p(s)  at  s  =  0} 

Then  a  set  of  compensators  c(s)  exists  that  will  simultaneously  stabilize 
the  feedback  loop  H  (s)  (without  RHP  pole-zero  cancellations  between 

Vl 

p(s)  and  c(s)),  and  meet  the  initial  condition  constraints 

{V2(0'^) .  ,v^(0+)}  =  iff: 

1.  r{t(s)}  +  r{p(s)}  £  M  +  1. 

and 


2.  (d.  (s),d  (s))  =  1. 

p 

Where  y^  =  0  for  i  <  M,  and  y^^  /  0. 


Proof 

The  lemmas  and  properties  derived  earlier  were  dependent  on  the 
assumption  that  r{p(s)}  £  0,  Then,  Theorem  3.1  and  Figure  4  lead  to 
Equation  3.12,  which  defines  the  complete  set  of  compensators  that 
meet  both  feedback  loop  stability  and  the  n-th  derivative  initial 
response  constraint  are  those  which  admit  stable  w(s)  and  z^(s).  Next, 
Property  3.10  states  that  both  w(s)  and  z^{s)  will  be  stable  iff  9^,(5), 
k^{s),  and  x^{s)  admit  stable  solutions,  as  described  per  Equations  3.33, 
3.34,  and  3.35.  Then,  Lemma  3.1  and  Lemma  3.2  give  the  necessary  and 
sufficient  conditions  to  stabilize  these  three  quantities  to  be 

(d^(s),dp(s))  =  1 


<3.119> 
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and 

r{t(s)}  +  r{p(s)}  <  M  +  1  <3.120> 

where  =  0  for  i  <  M,  and  f  0.  Thus,  <3.120>  leads  directly  to 
condition  1  in  the  theorem,  while  <3.119>  is  dependent  on  the  derivative 
constraint  considered.  Finally,  Lemma  3.4  states  that  if  Equations  3.100 
or  3.101  are  satisfied,  then  the  n-th  derivative  solution  also  stabilizes 
all  lower  derivatives.  Therefore,  to  get  all  constraints,  let  n  =  j 
for  the  constraints  {yQ,...,yj}  leading  to  the  assumptions  numbered 
1  and  2  in  Theorem  3.2,  as  well  as  d^(s)  =  d^(s)  in  Equation  3.119,  so 

that  (d^(s),d  (s))  =  1  is  required. 
z.  p 
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To  reiterate,  as  the  proof  for  Lemma  3.4  showed,  <3,119>  and  <3.120> 
may  not  be  necessary  assumptions,  but  were  shown  to  be  sufficient.  Also, 
obviously  these  pose  very  little  restriction  anyway,  so  that  this  point 
will  not  be  labored  further.  Now  the  paramterization  will  be  given  for 
the  compensators  that  achieve  these  constraints. 


Theorem  3.3  -  Initial  Response  Parameterization  (IRP) 

Given  that  the  assumptions  and  conditions  posed  in  the  Initial 
Response  Theorem  are  met  for  the  constriants  {v^CO'*’) , . .  .v^(0''’) }  = 
{yQ«--->yj}-  Then,  let  the  plant  be  comprimely  represented  by 
nn(s) 


p(s) 


so  that  stable  Up(s)  and  Vp(s)  exist  such  that 


u„(s)n„(s)  +  v„(s)d„(s)  =  1 


Also,  define 


where  n.(s)  and  d^(s)  are  coprime.  Next,  define 
j  j 


where  n^(s)  and  d^(s)  are  coprime,  and  let 


where  n.(s)  and  d.(s)  are  coprime.  Also,  define  stable  u  (s)  and 

V  (s)  such  that 
j 

Us(s)n  (s)d  (s)n.(s)  +  v  (s)d  (s)  =  1 
j  j  ^  j  j  j 


Then,  the  complete  set  of  compensators  that  stabilize  the  feedback  loop 
for  the  system  in  Figure  2  (without  pole-zero  cancellations  between 
p(s)  and  c(s)),  and  meet  all  j  +  1  initial  response  constraint  is  given 
by 


-w(s)d  (s)  +  u  (s) 

“  w{s)np(s)  +  Vp(s)" 

where  w(s)  is  stable,  and  given  by 


w(s)  =  {Up(s)ng(s)n^(s)  -  n^ 


u  (s)  +  d  (s)e(s) 


where  e(s)  is  an  arbitrarily  stable  function  such  that  w(s)np(s)  +  Vp(s) 
is  not  identically  zero. 


Proof 

Theorom  3.2  states  the  assumptions  and  resulting  N  &  S  conditions 
for  a  parameterization  to  exist.  Also,  Lemma  3.3  gives  that  complete 
parameterization  for  w(s)  leading  to  the  n-th  derivative  constraint  to 
be  met.  If  we  let  n  =  j  where  j  is  the  maximum  order  derivative  initial 
constraint,  then  the  assumptions  in  the  IRT  guarantee  that  for  n  =  j  in 
Equation  3.82,  all  constraints  will  simultaneously  be  met.  Also,  the 
proof  for  Lemma  3.3  showed  that  at  least  one  e(s)  exists  so  that 
w(s)np(s)  +  Vp(s)  is  not  zero,  in  which  case  the  resulting  compensator 
would  be  undefined,  even  though  z„(s)  and  w(s)  may  be  stable  by  Equation 
3.12.  Finally,  c(s)  was  defined  in  Equation  3.8  exactly  as  stated  in 
Theorom  3.3. 
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3.8  Summary 

The  IRP  can  be  used  to  find  the  complete  set  of  compensators 
that  will  meet  any  given  initial  response  conditions  that  are  within 
the  bounds  of  the  constraints  and  assumptions  given  in  the  IRT.  This 
then  gives  the  designer  a  tool  that  may  help  him  to  shape  the  transient 
response  of  a  system  as  he  desires.  For  example,  designing  for  large 
positive  initial  values  on  the  initial  response  derivatives  should  lead 
to  a  "faster"  response  than  designing  for  smaller  initial  values  on 
the  same  derivatives. 

I 


Theoretically,  these  theorems  could  be  used  to  design  precisely 
an  entire  system  response.  This  can  be  seen  by  representing  the 
frequency  response  of  the  entire  system  in  its  power  series  form,  as 
follows, 

H(s)  =  a  s'^  +  a,s"^  +  a«s~^  +  a,s"^  <3.121> 

0  Li 

Of  course,  this  transforms  to  the  time  domain  as 
h(t)  =  {ag  +  a^t  +  ^  t^  +  ^  t^  +...}U(t) 

CO  3  , 

=E  TT  t’u(t)  <3.122> 

i=0 

It  is  easy  to  see  that  the  initial  i-th  derivative  is  equal  to  a^. ,i.e., 

h^(0+)  =  a.  <3.123> 

But,  the  IRP  allows  one  to  specify  any  and  all  h’(0'^)  whenever  the 
bounds  of  the  IRT  are  not  crossed.  Thus,  for  any  i,  a^  can  be  specified, 
and  looking  back  to  <3.121>,  one  can  see  that  any  or  all  (theoroetically) 
of  the  power  series  for  the  transfer  function  can  be  specified! 


CHAPTER  IV 


EXAMPLES 

4.1  Introduction 

In  the  previous  chapter,  the  IRT  and  IRP  were  derived.  Now, 
several  examples  will  be  given  which  explore  the  usefulness  of  these 
theorems . 


4.2  Introduction  Example 

In  Section  1.2,  an  example  problem  was  given.  For  that  example. 
Equation  1.2  and  1.8  were  given  to  describe  the  set  of  compensators  that 
meet  the  desired  constraints,  but  the  equations  were  not  derived.  Thus, 
the  IRT  and  IRP  will  now  be  used  to  solve  this  problem.  Figure  1 
pictured  the  system  in  question,  which  is  given  again  in  Figure  5.  For 
the  system,  the  constraints  given  were  that  the  feedback  loop  be  stable, 
and  V2(0'’‘)  equal  3(=  y^) . 

For  this  system,  r{p(s)}  =  1,  and  the  input  has  fewer  than  three 
poles  at  s  =  0  (it  actually  has  none).  Therefore,  the  assumptions 
necessary  for  the  IRT  are  met.  Then,  since  yp  =  3,  M  is  equal  to  zero. 
Thus,  the 


V2(t) 


Figure  5.  Example  From  Introduction 


first  condition  given  in  the  IRT  is  met,  i.e.. 
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r{t{s)}  +  r{p(s)}  _<  M  +  1  ? 

0+l  =  l<0  +  l  =  l  / 

Also,  since  j  =  0,  Equation  2.39  leads  to 


t(,(s)  = 

This  can  be  coprimely  represented  by 
2 


to(s)  = 


(s+l) 


7 


1 


(s+l)' 


n.(s) 


The  plant  can  be  coprimely  represented  by 


1 


p(s)  =  ^ 


s-1 

s+r 


n„(s) 


Obviously,  (d^(s),d  (s))  =  1,  so  the  IRT  dictates  that  a  set  of 
tQ  P 

compsensators  exists  that  meets  the  desired  constraints. 


<4.1> 


<4.2> 


Now,  moving  on  to  the  IRP,  there  must  exist  stable  Up(s)  and  Vp(s) 
such  that 

Up(s)np(s)  +  Vp(s)dp(s)  =  1 

or 

Obvious  solutions  to  <4.3>  are 

Up(s)  =  2 

and  <4.4> 


v„(s)  =  1 


where  w(s)  will  be  described  momentarily.  Of  course,  <4.9>  is  exactly 


equal  to  <1.2>.  The  IRP  also  describes  w(s)  by 


djs) 


w(s)  =  {u  (s)n  (s)n.(s)  -  n  (s)  .  +  d  (s)e(s) 


0  "0 


where  e(s)  is  arbitrarily  stable  such  that  w(s)np(s)  +  Vp(s)  is  not 
identically  zero.  Substituting  in  values  for  the  current  functions 
leads  to 


w(s)  =  {2.1.—^  -  3s  ^}1  +  ^  e(s) 
(s+1)^  ^  ^  ^  ^ 

which  simplifies  to 


w(s)  =  —  —  I  +  ^  e(s) 
which  is  identical  to  <1.8>. 


4.3  dc  Motor  Initial  Response  Problem 


As  a  second  example,  consider  for  the  plant  a  field-controlled 
dc  motor.  Reference  (6)  gives  the  transfer  function  for  this  type  of 
motor  as 


"  V7  "  fen+i 


where  the  parameters  are  defined  as  follows, 
w  =  velocity  -  rad/sec 
v^  =  input  voltage  -  volts 
K  =  motor  constant  -  Ib-ft/field  amperes 
R  =  field  resistance  -  ohms 
F  =  total  viscous  friction  -  Ib-ft/rad/sec 


Intuition  dictates  that  controlling  even  just  the  value  of  w  (0  ) 
should  give  some  control  over  the  transient  response.  Thus,  this  theory 
will  be  tested  by  requiring  {w(0'’‘) ,w  {0+) ,w  (O'*') }  =  {0,0, A}  where  A  is 
a  real  number.  One  would  expect  large  values  of  A  to  lead  to  faster 
responses,  with  more  overshoot  than  small  values  of  A  generally. 

To  that  end,  we  apply  the  IRP  with  j  =  2,  and  assume  that  neither 
T^  or  T^  are  equal  to  zero  or  one.  Then,  a  coprime  representation  for 
p(s)  is 


_  (s+l)2 

'■  (i+sT^Hi+sT^)  - 


(s+l)‘ 


which  leads  to 


b,s+b- 

u„(s) 


<4.14> 


<4.15> 


Vp(s) 


b3S+b4 


where 


<4.16> 


^3  =  rf- 

'f  m 


<4.17> 


3-b3(T,.Tj 


Also,  a  coprime  representation  for  t2(s)  is 


Vs)  =  s'*  I  = 


n.(s) 

tg 


^  j - 1 - 

(s.l)3  ‘2 

and  r2{s)  is  coprimely  represented  by 


2  p  . , 1  As 


Equations  4.14  and  4.18  lead  to 


n„(s) 


‘2^^^  ■  dtlsT  1 

^  (3.1) 


=  (s+1)  ^  K_ 


(s+1) 


which  can  be  coprimely  represented  as 


jL  P;,(s) 

fO  -  RF  _  ^2 

®2^^^  "  1  "  dTsT 

TsTiy  ^2 


n-(s) 
_2 
d  (s) 
^2 


The  IVP  also  dictates  that  there  are  stable  u.{s)  and  v^(s)  such  that 


u  (s)n  (s)d  (s)n.(s)  +  v  (s)d  (s)  =  1 
^2  ^2  ^  ^2  ^2  °2 


and  substituting 


s''(l+sTj(l+sT„) 


Solutions 

for  u.(s) 
^2 

and  V  (s)  are 
S2 

c,s+c. 

u  (s) 
^2 

_  1  L 

S+1 

and 

4 

c-s  +c, 

3  2 

,S  +CcS  +CcS+C 

c  =  i— RL-T.y.y  <4  2 

^2  KV  T  T 
RT  Vm 

Note  that  c^-c^  can  easily  be  found,  but  will  not  be  written  here  since 
V  (s)  is  not  needed  to  solve  for  w{s).  Now,  the  entire  solution  set 

for  w(s)  is  solved  as  follows. 

w(s)  =  {u„(s)n  (s)n. (s)  -  r„(s)d^{s) }u,(s)  +  d,(s)e(s) 


d.s  +d„s  +d-s  +d.s  +dj-s  , 

^ — - — ^-5^ — ^  m 

(s+1)^  ^  ^ 


1  RF 

j  =  {b,-A  ;  ^2  through  dg  =  constants 


and  e(s)  is  arbitrarily  stable.  Thus,  the  set  of  stabilizing 

A 

compensators  that  cause  w^(0'*')  =  A  are  given  by 


c(s)  = 


-w(s)dp(s)+Up(s) 

w(s)np(s)+Vp(s) 


(l+sT,)(l+sT  )  b,s+bp 
-w(s)  - f 

s+U^  ^  ‘ 


bjS+b^ 


<4.28> 


(s+l)‘ 


where  {b^.bg.b^jb^}  are  given  by  <4.17>,  and  w(s)  is  given  by  <4.26> 
and  <4.27>.  Note  that  since  e(s)  is  arbitrarily  stable,  but  is  multiplied 
by  l/(s+l),  it  can  always  be  used  to  change  any  of  the  constants 
{d2,d2,d^,dg}  arbitrarily,  but  can  not  modify  d^.  Thus,  it  is  important 
only  to  know  dj. 

One  obvious  choice  for  e(s)  is  to  set  it  so  that 


d,(s+l)'’ 

w(s)  =  - ^  =  d, 

(s+1)^  ’ 


<4.29> 


For  example  purposes,  such  a  compensator  was  chosen  to  demonstrate  how 
varying  A  affects  the  transient  response  of  the  system.  Invoking  Equation 
3.10  leads  to  (where  h(s)  represents  the  feedback  loop) 
t(s)h(s)  =  t(s){-w(s)np(s)dp(s)  +  Up(s)  n^{s)} 


(l+sT^){l+sTj  (b,s+bo)  Sf 


SH 


which  after  some  manipulation  reduces  to 

^  {bj+b2-dj(T^  +  T|^)}s  +  ^  (b2-d^) 

t(s)h(s)  - - - -  <4.30> 

s{s+l)^ 

A  partial  fractions  expansion  eventually  helps  lead  to  the  time  response 
L~^{{t(s)h(s)}  =  th(t) 

(A-ep)  p  (e,-ep-A)  . 

=  {e2+e  ^(-e2-e2t  +  ^  - r)}U(s)  <4.31> 

where 

®1  ^  W  <4.32> 

and 

®2  "  W  (b2-d^)  <4.33> 

Now,  assume  that  the  input  v  =  100  volts,  and  the  motor  has  the 
following  constants. 

K  =  60  Ib-ft/amp 

R  =  50  ohm 

F  =  2  Ib-ft/rad/sec  4.34 

=  .4  henry/ohm 

T  =  8.75  rad-sec 
m 

These  values  were  used  with  several  values  for  A,  and  the  resulting 
time  responses  were  plotted  in  Figure  7  and  Figure  8.  Note  that  as 
expected,  even  though  only  w‘^(0''’)(=A)  was  directly  controlled,  increasing 


this  produced  faster  responses  and  eventually  increasing  overshoot. 


Response  for  A  =  10 


Figure  7.  Time  Response  Curves 


4.4  Combined  Constraints  Problem 


The  previous  problem  dealt  with  controlling  only  an  initial  value. 
Though  the  design  technique  yields  feedback  loop  stability,  there  were 
no  constraints  on  the  asymptotic  value  of  the  motor  velocity.  Thus,  as 
Figure  7  and  Figure  8  show,  using  the  previous  technique  alone  can  lead 
to  a  wide  variety  of  final  values.  Of  course,  a  common  practical  problem 
could  arise  in  which  both  the  initial  response,  and  the  asymptotic 
response  must  be  simultaneously  controlled.  Although  a  general 
parameterization  has  not  yet  been  derived,  specific  problems  of  this 
type  can  be  treated  now  by  invoking  both  the  IVP,  and  the  Tracking  Theorem 
given  in  (10)  and  Appendix  B.  Compensators  which  meet  both  theorems 
can  be  found,  and  then  analyzed  to  find  the  common  parameterization. 


As  an  example,  consider  the  same  motor  in  the  previous  section, 
with  the  same  constants,  and  the  same  constraint  w  (0  )  =  A.  Thus, 
<4.26>  describes  w(s),  and  <4.28>  describes  c(s).  Now,  assume  that  the 
velocity  of  the  motor  shaft  is  to  follow  the  input  voltage  v  so  that 
w(®)  =  V.  This  part  of  the  problem  can  be  solved  by  the  Tracking 
Theorem  iff  the  necessary  and  sufficient  condition  that  np(s)  and  d^(s) 
are  coprime  is  met.  Of  course,  this  would  not  necessarily  guarantee 
that  both  the  IRP  and  the  Tracking  Theorem  have  any  common  solutions. 
Checking  this  conditions,  t(s)  was  given  in  <4.12>,  and  can  be  coprimely 
represented  by 


Equation  4.14  gives  n  (s),  which  is  obviously  coprime  with  d^(s)  above, 
so  the  condition  is  met.  Thus,  according  to  the  Tracking  Theorem, 
there  exist  stable 

which  leads  to 


Unt(s) 


K 

RF 


•pt 


(s+1)' 


r  ^  i+r '  ‘ 


Solutions  for  <4.37>  are  easily  found  to  be 


and 


s+2 

s+1 


Also 


,,,  ,  v=’ 

"  -3^-- 


(s-H)‘ 

(l+sT^)(l+sTj 


s(s+l) 

which  can  be  coprimely  represented  by  a'(s) 


<4.37> 


<4.38> 


<4.39> 


n,,(s)/d^,(s),  where 

a  d 


na.(s)  = 


(l+sT^)(l+sT^) 


{s+l)‘ 


and 


<4.40> 


Now,  the  Tracking  Theorem  gives  the  solution  to  be 


{-w'(s)d  (s)  +  u  (s)} 

"  {wTslnJs)  +  Vp(s)} 

where 

w'(s)  =  -‘Jpt(s)Vp(s)  +  e'(s)d^.(s)  <4.42> 

with  e'(s)  an  arbitrarily  stable  function.  Therefore,  substituting 
into  <4.42>  gives 


w'(s) 


RF 

K  s+1 


+  e'(s) 


s 

s+1 


<4.43> 


Comparing  <4.28>  and  <4.41>,  one  can  conclude  that  simultaneous 
solutions  to  both  the  initial  value  problem  and  the  tracking  problem  are 
those  where  w(s)  =  w'(s).  Thus,  equating  <4.26>  and  <4,43>  leads  to 


d,s^+d,s^+d,s^+d.s^+d^s  , 


(s*l)® 

K  s+1  ®  s+1 

which  can  be  rearranged  to 


s+1 


e(s)  *  5+re'(s)  = 


d,s^+d,s^+d,s^+d„s^+dcS 
1  2  o  r  b 


(s+1)' 


RF 


<4.44> 


K  s+1 

The  above  equation  can  now  be  solved  be  differential  equation  methods, 
as  follows.  First,  a  particular  solution  can  be  found  by  recognizing 
the  coprimeness  of  the  functions  multiplying  e(s)  and  e'(s).  Proceeding 
with  this,  we  find  that 


s 


The  compensator  can  now  be  represented  by  <4.41>  and  <4.42>  or  by  <4.28> 
and  <4.26>.  Using  <4.42>,  we  find  that  w'(s)  is  given  by 


w‘(s)  =  w(s)  =  -Up^(s)Vp(s)  +  e'{s)  dg,(s) 

RF  '*’^3^ 

’  ■  ^  ="i  ”  (sn)= 

^RF  ^*^*4,^,,)  1.  s  ^4, 

K  5+1  '  5+1'  5+1  ■ 

As  in  the  previous  example,  a  particular  k(s)  has  been  chosen  for 

example  purposes.  In  this  case,  k(s)  can  be  chosen  so  that  the  term 

inside  the  {  }  in  <4.51>  is  equal  to  dj  +  (RF/K)b2.  Using  this  k(s) 

./.i  _  RF  ,  RF  ,  ,  s 

■  ■  K  s+1  ^^1  K  ^3^  s+1 


He  Rf"  h 
-  K  ^4 


Now,  invoking  <3.10>  leads  to  (where  h{s)  represents  the  feedback  loop) 


V  W  (l+sTj(l+sT) 

t(s)h(s)  - 1. - ^ 

^  (s+1)^  (s+1)^ 


.  W  4 


which  reduces  to 


3  2 

As  +q,s  +q-s+v 


and 


■>2  '  W  *  vb,TfT„ 


<4.54> 


>’3  '  W  *  ''>=4<Tf 

Invoking  the  Final  Value  Theorem  on  <4.53>,  we  find  that  asymptotical ly 
w(t)  is  equal  to  v.  Also,  the  Initial  Value  Theorem  shows  that 
w  (0  )  =  A.  Thus,  we  do  in  fact  have  an  adequate  compensator.  Now,  a 
partial  fractions  expansion  leads  to  the  time  domain  response 

L“^  {t(s)h(s)}  =  th(t) 

2  3 

=  {v+e’^(-v-vt+(A-v)  +  (q2-sA-v)  |- 

.4 

+  (A-v+q3-q2)  ^  }U(t)  <4.55> 

The  constants  in  Equation  4.34  were  used  for  Equation  4.55.  The  resulting 
response  curves  were  again  plotted  for  several  values  fo  A  in  Figure  9 
and  Figure  10.  Note  that  the  results  are  again  as  expected.  Increasing 
A  generally  increases  the  speed  of  the  response,  and  the  amount  of 
overshoot,  while  all  of  the  curves  have  the  steady  state  value  w  =  100 
rad/sec. 
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a.  Response  for  A  =  -100  and  v  =  100 


b.  Response  for  A  =  0  and  v  =  100 


Figure  9.  Time  Response  Curves 


CHAPTER  V 


RECOMMENDATIONS 

The  work  in  this  thesis  leaves  the  doors  open  for  the  study  of. 
several  related  topics.  The  following  list  contains  recommendations 
of  a  few  such  problems  which  could  be  researched. 

1.  The  problem  considered  here  only  utilizes  unity  gain  in  the 
feedback  loop.  An  obvious  extension  is  to  consider  the  more  general 
case  with  non-unity  feedback  gain. 

2.  To  be  truly  useful  for  "real  world"  applications,  this  single¬ 
variate  case  should  be  extended  to  multivariate.  Similar  extensions 
have  previously  been  handled  for  other  parameter! zations,  for  example 
in  (1,  4,  7,  14,  16,  17,  18). 

3.  A  third  extension  would  be  to  develop  a  simultaneous  Initial 
Response-Tracking  parameterization,  as  well  as  other  similar  simultaneous 
parameter! zations.  These  would  be  based  on  the  IRT  and  IRP,  and 
parameter! zations  such  as  those  in  (10,  11). 

4.  The  examples  in  the  previous  chapter  were  designed  to  illustrate 
the  usefulness  of  specifying  initial  response  values.  This  technique 
should  be  explored  and  experimented  with  further  so  that  perhaps  a 
general  design  method  would  be  derived.  The  power  of  specifying  multiple 
initial  response  values  particularly  needs  further  exploration. 

5.  At  the  end  of  Chapter  III,  it  was  illustrated  that  specifying 
initial  response  values  is  equivalent  to  specifying  the  coefficients 
of  the  power  series  of  the  transfer  function.  To  be  completely 


specified,  of  course,  an  infinite  number  of  coefficients  must  be 
specified.  However,  the  concept  of  approximately  specifying  the  system 
transfer  function  by  only  specifying  a  finite  number  of  coefficients 
could  be  pursued.  This  could  lead  to  the  ability  to  design  for  an 
entire  specific  response  curve  simply  by  applying  the  IRT  and  the  IRP. 

6.  The  IRT  contains  two  assumptions  based  on  Lemma  3.4  which 
guarantee  that  the  parameterization  for  a  particular  initial  response 
derivative  will  also  meet  the  lower  order  derivative  constraints.  The 
assumptions  are  not  very  limiting,  but  it  is  possible  that  they  are  not 
required  at  all,  or  that  less  stringent  assumptions  could  be  used. 

This  problem  could  be  further  researched. 

7.  The  results  in  this  thesis  depend  upon  the  system  model,  and 
no  consideration  has  been  given  for  the  robustness  necessary  usually 
in  the  real  world  due  to  approximations,  drift,  etc.  This  robustness 
topic  should  be  looked  into  further,  as  it  has  been  for  other 
parameterizations ,  as  in  (11). 
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Several  properties  were  given  in  Section  2.3  that  will  be  proved 
in  this  appendix.  The  properties  relate  to  the  transformation  described 


in  Equation  2.17,  which  is  repeated  below. 


H(s)  =  T{H(s)}  =  H(i) 
T'^{H{s)}  =  H(^) 


Additive  Property 

H(s)  +  F(s)  =  T{H(s)  +  F(s)} 


<A.1> 


<A.2> 


Proof 


Represent  any  two  transfer  functions,  H(s)  and  F(s)  as  follows, 


H(s)  = 


(ajS+bj)...(a^s+b^) 


'1  1' 


m  m' 


(ejS+fj)...(e  s+f  ) 


gjS+jj)...(g^s+j 


<A.3> 


Of  course,{a^. ,b^,c^,d^. ,e^. ,f^,g^,j\}  are  complex  numbers,  and  {n,m,p,k} 
are  positive  integers.  Then,  according  to  <A.1>, 


H(s)  =  H(i) 


_  (r  ^  b^)...(/+  b^)  ^ 

Cl  ’  ^n 

_  (ai+bjs)...(a^+bj^h"’"'^ 

=  (ci+disj...(c^+d^) 


<A.4> 


°  — : — 9ir~  ■ 

(j  +  +  j|^) 


'yj  ■ 

Adding  <A.4>  and  <A.5>  yields 


(ej  +  fil)...(ep  +  fpS)i 

(gi  +  jiS.)...(gk  +  jbl) 


H(s)  +  F(s)  = 


(aj  +  bj^)...(a^  +  b^)s^ 

(c,  +dis)...(c„  +  d^) 

1  1—'  '  m  m— ' 


(ej  +  fil)...(e  +  Oh' 


.gi  +  ■'■  ^'k-^ 


Conversely,  by  <A.3>, 


(a,s  +  b,)...(a  s  +  b  ) 
(djS  +  fi)...{e  s  +  f  ) 


(g^s  +  jj)...(g|^s  + 
Then,  transforming  the  sum  leads  to 


T{H(s)  +  F(s)}  = 


jj...(r^+  jjs' 


(a,  +  b,s)...(a„  +  b  s)s"-” 
1  1—  m  m— 


(ei  +  fjl)...{ep  +  fpS)i 

(gi  +  jii)...(gk  +  j[l) 


<A.8> 


Equating  <A.6>  and  <A.8>  leads  directly  to 


H(s)  +  F(s)  =  T{H(s)  +  F(s)} 


<A.9> 


Multiplicative  Propertv 


H(s)F(s)  =  T{H(s)F(s)} 


<A.10> 


Proof 


This  proof  proceeds  exactly  like  the  proof  for  the  Additive  Property, 
except  that  the  product  of  H(s)  and  F(s)  is  considered  rather  than  the 
sum.  Hence,  this  proof  will  not  be  reiterated. 


Inverse  Transform  Propertv 


H(j)  =  H(s) 


<A.11> 


Proof 


Again,  represent  H(s)  per  <A.3>,  and  the  resulting  H{s)  per  <A.4>. 


b,  b 

-  1  (a,  +  T^)...(a^  +  -^) 


,  (»!=  ^  b,)...(a„s  t  b„)s 

(a^s  +  bj)...(a^s  +  b^) 
(CjS  +  dj)...(c^s  +  d^) 

=  H(s). 

EOP 


Transformation  Stability  Property 

If  H(s)  is  stable  (or  unstable)  in  the  s-domain,  then  so  is  H(£) 
in  the  s-domain,  and  visa  versa. 


Proof 


Once  more  define  H(s)  generally  by  Equation  A. 3  so  that 


_  (a,s  tb,)...(v^b„) 
(CjS  t  dj)...(c^s  1-  ij 


and  the  resulting  H(s)  by  <A.4>,  thus 


(a,  +  b,s)...(a  +  b^)s 

=  (cSd7)...(c>d5 


m-n 


Then,  first  consider  only  stable  H(s).  Under  this  condition,  n  j<  m, 
and  the  poles  are  all  in  the  strict  LHP.  Of  course,  to  meet  this  LHP 


requirement,  the  poles,  p^ ,  which  are  located  at 


must  be  such  that 

[p^ [  >0,  but  finite 


<A.12> 


<A.13> 


91 

and 

0.:  90°  <  0^.  <  180°,  or  -180°>  0^.  >  -90°  <A.l' 

To  meet  Equation  A. 13  requires  that  {c^.  ,d^.}  contains  no  elements  equal 
to  zero,  but  no  such  requirement  is  on  {a^. ,b^. }.  Therefore,  if  there 
are  A  elements  in  the  set  {b^},  1  <.  i  £  n,  equal  to  zero,  then 

r{H(£)}  =  m  -  (n-A)  -  (m-n)  =  A  ^  0  <A.l 

proving  that  H(£)  is  proper,  as  required  for  stability.  Also,  H(£) 
has  its  poles  located  at  p^ ,  where 

Pi  =  -  -  IPiliSl 

Since  lp^|  is  greater  than  zero  and  less  than  infinity,  so  is  |p. [. 
Also,  since  0^.  =  -0^,  then  by  <A.14>, 

0.  :  90°  <  0.  <  180°,  or  -180°  >  0.  >  -90°  <A.l' 

so  that  0^.  is  in  the  LHP.  Therefore,  H(£)  is  stable  if  H(s)  is. 

Consider  now  unstable  H(s).  First  consider  the  case  where  the 
transfer  function  is  net  proper,  so  that  n  >  m.  Analysis  of  <A.4> 
illustrates  that  this  results  in  (n-m)  poles  at  s_  =  0  for  H(£),  which 
makes  H(s^)  unstable.  Next,  consider  finite  poles  of  H(s)  in  the  closed 
RHP.  In  this  case 

0.  :  -90°  <  .j  <  90°  <A.ll 

which  by  Equation  A. 16  results  in  0^  in  this  same  unstable  region, 
since  0^  =  -0^..  Therefore,  unstable  H(s)  lead  to  unstable  H(£). 

Conversely,  the  inverse  transform  H(^)  is  equal  to  H(s)  by  the 


Reference  (10)  states  the  following  theorom.  The  proof  can  be 


found  in  that  reference. 

Tracking  Theorem 

Give  p(s)  there  exists  a  compensator  for  the  feedback  system  of 

Figure  2  which  stabilizes  the  feedback  loop  and  simultaneously  causes 

it  to  track  the  impulse  response  of  t(s)  if  and  only  if  n  (s)  and  d.(s) 

P  ^ 

are  coprime.  In  this  case  let  Up^(s)  and  Vp^(s)  be  stable  functions 
such  that 

Upt<s)np(s)  t  Vpp(s)dp(s)  =  1 

and  let  a‘(s)  =  n  ,(s)/d,,(s)  be  a  coprime  fractional  representation  of 

d  a 

3'(s)  =  d  (s)/d.(s).  Then  the  desired  set  of  compensators  take  the  form 
P  ^ 

{-w(s)d  (s)  +  u  (s)} 

j.(s)  -  {v^(s)np(s)  +  Vp(s)} 

where 

^  e(s)dg,(s) 

with  e(s)  an  arbitrarily  stable  function  such  that  w(s)np(s)  +  Vp(s) 
is  not  identically  zero. 
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